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Eliminating Error in the Chemical Abundance Scale for 
Extragalactic H II Regions 
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ABSTRACT 

In an attempt to remove the systematic errors which have plagued the calibration of 
the H II region abundance sequence, we have theoretically modeled the extragalactic 
Hll region sequence. We then used the theoretical spectra so generated in a double 
blind experiment to recover the chemical abundances using both the classical electron 
temperature + ionization correction factor technique, and the technique which de- 
pends on the use of strong emission lines (SELs) in the nebular spectrum to estimate 
the abundance of oxygen. We find a number of systematic trends, and we provide 
correction formulae which should remove systematic errors in the electron tempera- 
ture + ionization correction factor technique. We also provide a critical evaluation of 
the various semi-empirical SEL techniques. Finally, we offer a scheme which should 
help to eliminate systematic errors in the SEL-derived chemical abundance scale for 
extragalactic Hn regions. 

Key words: galaxies: abundances, ISM: abundances, Hll regions, methods: data 
analysis, techniques: spectroscopic 



1 INTRODUCTION 

To understand how galaxies in the early universe evolved 
into those that we see locally requires an understand- 
ing of both chemical and star formation history of galax- 
ies over cosmic time. Indeed, observational data com- 
piling hundreds of thousands galaxies at different red- 
shifts are suggesting that there actually is a deep phys- 
ical connection b etween metallicities, masses and star- 
formation rate (e.g. Tremonti et al. 2004; Brinchman net al.l 



l2004lKewlev fc Ellisonll2008l ; lLara- Lopez et al.ll2010al lbh, al- 
though these results are mainly based on trends of the chem- 
ical abundances and not on their absolute metallicities. 

Theoretical simulations can now p redict the evolution 
of chemical abundances in g alaxies (e.g. lMolla fc Diazn2005l; 
iDave fc Oppenheimerl 120071 ). However, it is difficult to put 
these predictions to the test because the absolute calibration 
of the chemical abundance scale from H II regions is currently 
discrepant by more than a factor o f two. See, for example 
the work of lKewlev fc EllisorJ (2008) where abundances were 
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calculated for 45,000 galaxies fr om the Sloan Digital Sky 
Survey (SDSS) (|York et al.ll2000l : lAbazaiian et al. .1120031 ) us- 
ing 10 different calibrations given in the literature. Although 
all methods show some correlation with abundance, there is 
typically a scatter ranging over a factor of 2-3 in the absolute 
oxygen abundance derive d from the various techniques. Th e 
same result was found bv lLopez-Sanchez fc EstebarJ IJ2010l ). 
who also considered in their analysis the metallicities de- 
rived using a direct determination of the electron tempera- 
ture of the ionized gas. Without a proper estimation of the 
absolute chemical abundances in galaxies, the comparison 
between observational data and theoretical models is not 
entirely satisfactory. 

Thanks to its strong emission lines, the easiest element 
to measure in Hll region emission spectra is oxygen. This 
is fortunate, since O is an a-process element made directly 
in short-lived massive stars. It comprises about 50% of all 
the heavy elements by mass throughout the universe, and is 
therefore representative of all the heavy elements. By con- 
trast, Fe is made in lower mass stars, and is not released 
promptly into the interstellar medium of galaxies. The O/H 
abundance in Hll regions -typically expressed in terms of 
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12+log(0/H), where O/H is the ratio of the number of oxy- 
gen to hydrogen atomic density . This has been determined 
in three semi-empirical ways, as well as through direct the- 
oretical photoionization modeling: 

1. By direct measurement of the electron tempera- 
ture, from line ratios such as the [O in] A4363/A5007 
(see lOsterbrock fc Ferlandl [2006) , followed by an anal- 
ysis of the various ionization fractions in the zones 
of the Hll region which produce optical emission 
lines. We dub this the T c method. This has been 
very extensive ly applied by observers over the past 
30 years ( e.g. Peimbert fc Costerd 1963; StasiriskalL 1978i; 
Ailed Il984l: iDfaz et all 1 19871: llzotov. Thuan. fc Lipoyetskil 



1994 lEsteban fc Peimbertl 1 19951 : IVflchez fc Estebanl 
Stasihskal 120051 : iBresolin et alj 120051: 



Izotov et a] 



199 



200 



Esteban et al.ll2004 lLopez-Sanchez fc Estebanl l2009h 

2. By direct measurement of recombination lines of ions of 
heavy elements such as O II and, in some cases, O I and com- 
paring these with the recombina t ion lines of hydrogen (e.g . 
Esteban et ail Il998l. 12004 I2009J; iGarci'a-Roias et al.1 1 2004 : 



Garcfa-Roi as & Est ebanl I2007J : lLopez-Sanchez et all 12007 : 
Peimbert et alj 120071 ). We refer to this technique as the RL 
method, although it is generally not useful in the determi- 
nation of the cosmic abundance scale, since the recombina- 
tion lines of heavy elements are generally too weak (~ 10 4-6 
times fainter than H/3) to measure. 

3. By measurement of the strong emission lines in the H n 
region spectrum calibrated using 

a. ph otoionizatio n mo de ling (e.g. | Paeel et al.l 



19791; IMcGaughl [l99ll; JKewlev fc Dopital 



Kobulnickv fc Kewlevl 120041 : iDopita et alT 



Dors fc Copettill2006l ) 



b. Hn regions for 
known from the T e 



Denicolo. Terlevich fc Terlevichl 12002 



which the oxy gen abundance is 
method (e.s. Pilvuginl l2001allbl; 



2004 



Pettini & Pagel 



Pilvugin fc Thuanl l2o"o5l; IPerez-Montero fc Diaz 



2003 ; iPilvugin et al.ll2oia IPilvugin fc MattssonlbOllI ) 



We will call this the Strong Emission Line (SEL) 
method, and it is nowadays extensively used to determine 
metallicities in large galaxy surveys. 

Each of these methods has its particular problems, and 
range of application. More importantly, where all methods 
can be used, the derived O/H abundances often show sys- 
tematic disagreements between the various methods. These 
may amount to factors of two or more. In this paper, we seek 
to resolve some of the causes of these discrepancies by re- 
deriving the chemical abundances of model H n regions using 
the T e method and the standard SEL techniques. We use a 
double-blind approach, so that the analysis of the photoion- 
izaton models was conducted without any a priori knowl- 
edge of the chemical abundances and physical parameters 
used in the photoionization models. 

We will now briefly indicate how these alarming abun- 
dance discrepancies between the various methods may arise. 

1.1 The T e and RL methods 

At normal nebular temperatures (T e ~ 10 4 K), the auroral 
[O ill] A4363 line is only a few percent as strong as the neb- 
ular A5007 line, so that high quality spectra are required to 



give an adequate signal-to-noise ratio to measure the fainter 
line. When the S/N ratio is inadequate, the strength of the 
fainter line tends to be overestimated, leading to an underes- 
timate in the derived abundance. This is also true for other 
commonly-used temperature sensitive line ratios. 

A more insidious pro blem was first pointed out by 
iPeimbert fc Costerd (|1969i ). This still has not been resolved 
to everyone's satisfaction. Normal Hll regions are not ho- 
mogeneous but contain temperature gradients, dense inclu- 
sions in which collisional de-excitation of cooling lines leads 
to higher temperatures, as well as colliding supersonic flows 
in which both density and electron temperatures may be 
raised as a result of shocks. In all such regions of enhanced 
temperature, the [O ill] A4363/A5007 emission line flux ratio 
is raised by a factor which depends only on the tempera- 
ture. This would not be a problem except for the fact that 
the line emissivity is also raised in such regions by a fac- 
tor proportional to the square of the electron density, nj,. 
Thus, the temperature estimate provided by the forbidden 
line ratio is dominated by the line ratio characterizing these 
over-dense inclusions, rather than providing a measure of 
the electron temperature of the nebula as a whole. The elec- 
tron temperature is therefore systematically overestimated. 
Hence, the presence of temperature gradients or fluctuations 
in the ionized gas leads to the abundances bas ed on the 
T e method being sy stematically underesti mated IjPeimberti 
1 19671 ; IPeimbert et~aT1l2007l ; IStasihskall2002l . 120051 ) . 

Indeed, detailed studies that compared heavy ele- 
ments abundances derived using b oth the T e method 
and the RL method in Galactic lEsteban ct al. 2004; 
2005, 2006; Garcfa-Roias & Esteban 
2007J; iMesa-Delgado et al.ll2009bl; iMesa-Delgado fc Esteban 



20ld; Mesa-Delgado et ah! 12011 ^ and extragalactic 



dTsamis et alTT2003l; [Peimberd bo!)! lLopez-Sanchez et al.1 
120071; lEsteban et all WQ^j ] Hn regions find a very good 
agreement between both results when temperature fluctu- 
ations are considered. However, their origin and existence 
is still controversial because they are not well reproduced 
by s t andard ph otoionizatio n models dKingdon fc Ferlandl 
1995; Rodriguez & Garcfa-Roias 2010), and hence addi- 
tional mechanisms are proposed to explain the presence of 
temperature fluctuations in the ionized ga s (see reviews by 
lEstebanll200J IPeimbert fc Peimbertll2006l '). 



1.2 The SEL method 

T he strong em i ssion line (SEL) method was first proposed 
bv lPagel et al.l 1)19791 ). It relies upon the ratio of the sum of 
the strong forbidden oxygen emission lines with respect to 
H/3, the so-called R23 ratio; R23 = I { [0 11] AA3726,3729 + 
[Oin] AA4959,5007) / /(H/3). Since then many other such 
ratios have been pr oposed. Th e most widel y used of these are 
the ca librations of IMcGaughl |l99lj) and iKewlev fc Dopital 
(2002), based on detailed photoioniza tion mod e lling, a nd 
the empirical rela t ions provid ed by iPilvuginl (feOOlallbr); 
IPilvugin fc Thuanl (|2005l ) and IPilvugin. Vilchez fc Thuanl 
1 20101) . Both kinds of calibrations strive to improve the ac- 



curacy by making use of the [O ill] /[On] ratio as ionization 
parameter, which accounts for the large scatter found in the 
-R23 versus oxygen abund ance calibration, which is larger 
than observational errors (|Kobulnickv et al.lll999l ) 

Hence, the SEL method is ideally suited to global abun- 
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Figure 1. The complete model grid plotted on the IVeilleux fc O sterbr ockl (119871^ di agnostics. [Oiii]/H/3 vs. [Nil] /Ho; (panel a) and 
[Om]/H/3 vs. [Sii]/Ha (panel b). In the lower panels we show the lDopita et al.l d200(J) diagnostics [Om]/ [On] vs. [Nn] /[On] (panel c) 
and [Om]/H/3 vs. [Nil] /[Oil] (panel d). The models are identified by their input (gas + dust) 12 + log(0/H) (written in bold face and 
color coded) and their input ionization parameter, logg. Models are stepped by 0.25 in logg. 



dance determinations in distant galaxies, since only the 
strong lines are visible against the stellar continuum. The 
abundance scale for ensem bles of H H reg i ons in galaxies has 
recently been calibrated bv lDopita et al.l (|2006al ). The basic 
problem with all ratios of optical forbidden lines to hydro- 
gen recombination lines is that they are two-valued as a 
function of chemical abundance. This abundance ambiguity 
can only be raised by the simultaneous use of several such 
ratios. An additional problem is that not all SEL method 
give similar oxygen abundances: the SEL methods based 
on calibrations using photoionization models generally tend 
to overpredict (but n ot always, see IPerez-Montero et al.l 
l20ld ; iDors et aUl201ll ') the observed oxygen abundances de- 
rived using the T c method by 0.2-0.4 dex d Pcimbert et al.l 
20071; iBresolin et alj 120091 ; lLopez-Sanchez fc Estebanll2oioT 
Moustakas et alj|20ich . As we should expect, that is not 



the case of the SEL techniques that are based on cal- 
ibrations with objects for which the oxygen abundances 
are well known from the T c method. Generally speaking, 



the SEL method generally returns abundances which are 
higher than the T e method, but in better agreement with 
those derived from recombinatio n lines of heavy elements 
JLopez-Sanchez fc E stcban 20ld '). 



1.3 A double-blind test 

In this paper, we seek to eliminate some of the system- 
atic errors in abundance determinations by applying the T c 
method and the SEL method to theoretically-generated H n 
region spectra, rather than using observed spectra. This has 
the advantage of eliminating observational errors but, more 
importantly, it also provides a direct test of whether the vari- 
ous methods can recover the chemical abundances and phys- 
ical conditions which were fed into the theoretical models. In 
order to ensure objectivity, Lopez-Sanchez (who applied the 
T e and the SEL method) was provided a randomized selec- 
tion of theoretical H n region models generated by Dopita. 
The data provided consisted of a set of line fluxes (relative 
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to H/3) as delivered by the model. The models were iden- 
tified only by letters, and no information was given about 
the input physical parameters which were input to gener- 
ate the models. This is the classical double-blind technique 
often applied in medicine, but rarely in astrophysics. 

In Section 2, we provide details of the techniques used 
to generate the theoretical spectra, in Section 3 we apply the 
T e to the theoretical spectra and examine how well the input 
parameters have been recovered. Where possible, we provide 
correction formulae to apply. Section 4 presents the results 
provided by the most common SEL techniques. In this sec- 
tion, we also compare the oxygen abundances provided by 
the SEL methods with the derived T c and RL abundances 
in high-quality Galactic and extragalactic Hll spectra. Fi- 
nally, Section 5 gives the conclusions and provides our rec- 
ommended procedure needed to accurately derive the oxy- 
gen abundance in Hn regions. 



2 THE THEORETICAL MODELS 

We have generated a sequence of photoionization models us- 
ing the Mappings his code, an updated vers ion of the code 
originally described in ISutherland fc Dopital (119931'). T hese 
model s are very similar to those produced by iDopita et al.l 
IJ2000h . and subsequently used bv IKewlev fc Dopital (|2002r i. 



but the chemical abundances have been specifically modified 
so as to provide a better fit to the Sloan Digital Sky Sur- 
vey (SPSS) observe d sequence of extragalactic H II regions 
iJKewlev et al.ll2006r i. 

Our grid of photoionization models uses as stellar 
input the fluxes com puted by the Starburst 99 code 
(jLeitherer et al.lll999T ) appropriate for continuous star for- 
mation. The input spectrum and the assumed stellar IMF 
inp ut into the Starbu rst 99 code is the same as used 
by IDopita et al.l J2000h. Th e later Starburst 99 models 
l|Vazauez fc Leithererll2005l ). which incorporate a fully self- 
consistent mass-loss formulation were not used as these gen- 
erate a stellar radiation field that is somewhat too so ft rel- 
ative to these earlier models (see IDopita et al. 2006a, for a 
detailed discussion of the issue). To summarise, the most 
likely cause of this discrepancy is probably that the stellar 
wind is clumpy as a result of radiation pressure instabili- 
ties, so a multi-dimensional model atmosphere would need 
to be constructed. Effects of stell ar rotation can also m odify 
the emergent ionizing spectrum IjLevesque et al.ll2010l ). For 
consistency, we have interpolated the stellar fluxes to match 
the chemical abundances actually used in the model, allow- 
ing for the abundance sc ale shift caused by th e re-calibration 
of the solar metallicity (jGrevesse et al.ll2010n . 

The extragalactic Hll region sequence is clearly 
defined in the observational plane by the standard di- 



agnostic line ratio plots of Baldwin. Phillips fc Terlevichl 
j 19811 ) and IVeilleux fc Osterbrockl (|l987l 1. These plot 
the ratios [N n] A6584/Ha, [Sn] AA6717, 31/Ha and 
[Oi1A 6300/Hq against th e [O ill] A5007/H/3 emission line 
ratio. iKewlev et al.l (|2006r ) gave these diagnostics for some 
45,000 galaxies drawn from the SDSS. Our photoionisation 
models closely reproduce the observational sequence of 
extragalactic Hll regions, bearing in mind that the ioniza- 
tion parameter tends to become higher as the abundance 



decreases, thanks to the dec reasing importance o f the ram 
pressure of the stellar winds IJDopita et alj|2006ah . 

The solar metallicity for our photoionisation mod- 
els was defined by t he solar oxygen abundance given by 
iGrevesse et all l|2010h . 12+log(O/H)=8.69±0.05. However, 
the abundance of helium and carbon were assumed to vary 
according the formula given bv lDopita et al.l (12006a ) . For ni- 
trogen, we use a form similar to that given in IDopita et al.l 
l|2006bh : 

12 + log(N/H) = 7.6 + log[Z/Z e + (Z/Z Q ) 2 }. (1) 

Here the zero point of the nitrogen abundance has 
been adjusted to prove a b etter fit in the corresponding 
IVeilleux fc Osterbrockl |l_987) diagnostic diagram. Because 
nitrogen is an important coolant at high metallicity, this 
has the incidental effect of improving the fit of the models 
to the observed extragalactic H II region abundance sequence 
in the other two diagnostic diagrams. For all models we in- 
cl ude dust, and dust p hysics according to the description 
in lDopita et al.l l|2006al ). where the depletion factors of the 
various elements from the gaseous phase are also given. 

Our models form a grid in chemical abundance or metal- 
licity, defined by the quantity 12+log(0/H), and the ioniza- 
tion parameter, logg, the ratio of the number of photons 
passing per unit area and time to the number density of 
hydrogen atoms. The q parameter is related to the dimen- 
sionless ionization parameter U by U = q/c, where c is the 
speed of light. 

Our photoionization models are spherical, with a target 
ionization parame t er for the inner surface. This is unlike the 
IKewlev fc Dopital |2002') models, which were plane parallel. 
Our models are also isobaric with log(P/fe) = 5.5 cm -3 K, 
where k is the Boltzmann's Constant, so that the density 
varies throug hout the model, but is typically 10 — 30 cm -3 . 
Unlike the IKewlev fc Dopital (2002) models, our models al- 
low for photoelectric heating of the gas by dust grains, and 
include the effects of radiation pressure, which is becoming 
appreciable at the higher values of the ionization parameter. 



2.1 The Double-Blind Test Models 

From the theoretical model Hll region grid we have taken, 
for a fixed metallicity, every second model in log q, and ap- 
plied to these the standard observational techniques for de- 
riving abundances. In order that no bias was incurred in 
the analysis, the spectra were randomly selected from the 
input grid of models, and the raw spectra labelled sim- 
ply A,B,C....Z, AA..AG were passed by Dopita to Lopez- 
Sanchez for analysis. This analysis returned physical condi- 
tions, ionization parameters, and chemical abundances given 
by each of the semi-empirical techniques commonly used by 
observers, which included both the T c and the SEL meth- 
ods. We can then compare the results of this analysis to 
the known inputs into the original models to discover where 
there exist systematical errors, and to evaluate the magni- 
tude of the random errors that can be generated by these 
techniques. The advantage of this approach is that it totally 
eliminates errors generated by the observations, and pro- 
vides spectra of essentially infinite signal to noise for analy- 
sis. 

In an ideal world, the parameters returned by the abun- 
dance analysis would be the same as those input into the 
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Table 1. The physical and chemical properties of the photoionization models used in the double- blind abundance fitting, organized 
according to chemical abundance 12+log(0/H), the nominal ionization parameter input into the model, the ionization parameter at the 
inner surface of the ionized gas, </; n , and mean ionization parameter in the ionized region of the model, < ?(H) >. For each model we give 
the emission line weighted mean temperatures (in K) for the [Nil], [Oin], [So] and [Oil] lines, and the mean recombination temperature 
for ionized hydrogen, T(rec). 



Name 


12+log(0/H) 


log q 


9i n 


< 9(H) > 


< T [N II] > 


< T [0 m] > 


< T [s „, > 


< T [0 ii] > 


T(rec) 


X 


9.39 


8.00 


5.417E+07 


2.598E+07 


3363 


2987 


4277 


4284 


3439 


J 


9.39 


7.50 


2.510E+07 


1.171E+07 


3041 


3620 


4201 


4194 


3908 


11 


9.39 


7.00 


1.005E+07 


4.426E+06 


2881 


4282 


4352 


4357 


4307 


c 


9.39 


6.50 


3.393E+06 


1.376E+06 


2617 


4568 


4494 


4558 


4403 


G 


9.17 


7.75 


5.358E+07 


2.432E+07 


4478 


4619 


5544 


5487 


5004 


AA 


9.17 


7.25 


2.140E+07 


9.721E+06 


4111 


5347 


5633 


5586 


5511 


O 


9.17 


6.75 


7.393E+07 


3.252E+06 


3766 


5765 


5745 


5782 


5745 


1 


8.99 


8.00 


1.212E+08 


4.573E+07 


5988 


5585 


6780 


6688 


5928 


R 


8.99 


7.50 


4.210E+07 


1.895E+07 


6014 


6162 


6801 


6688 


6433 


AC 


8.99 


7.00 


1.494E+07 


6.827E+06 


5826 


6736 


6898 


6872 


6836 


Z 


8.99 


6.50 


4.883E+07 


2.119E+06 


5351 


6936 


6865 


6935 


6902 


N 


8.69 


7.75 


1.025E+08 


3.850E+07 


8604 


8069 


8782 


8713 


8294 


A 


8.69 


7.25 


3.265E+07 


1.460E+07 


8680 


8448 


8751 


8729 


8618 


AG 


8.69 


6.75 


1.076E+07 


4.837E+06 


8738 


8764 


8693 


8791 


8780 


U 


8.39 


8.00 


2.720E+08 


6.599E+07 


10490 


10492 


10173 


10610 


10520 


F 


8.39 


7.50 


7.609E+07 


2.951E+07 


10414 


10492 


10124 


10517 


10490 


W 


8.39 


7.00 


2.314E+07 


1.030E+07 


10290 


10538 


10000 


10416 


10420 


AF 


8.39 


6.50 


7.096E+06 


3.055E+06 


9960 


10282 


9663 


10088 


10040 


T 


8.17 


7.75 


1.720E+08 


4.959E+07 


11335 


12146 


10784 


11676 


11960 


S 


8.17 


7.25 


4.807E+07 


1.997E+07 


11207 


11875 


10784 


11501 


11630 


L 


8.17 


6.75 


1.410E+07 


6.185E+06 


10815 


11379 


10353 


11060 


11020 


Q 


7.99 


8.00 


3.853E+08 


6.985E+07 


11843 


13584 


11123 


12474 


13250 


AE 


7.99 


7.50 


1.005E+08 


3.513E+07 


11792 


13180 


11140 


12366 


12820 


AB 


7.99 


7.00 


2.787E+07 


1.207E+07 


11483 


12420 


10897 


11912 


11990 


E 


7.99 


6.50 


7.929E+06 


3.320E+06 


10793 


11434 


10284 


11077 


10980 


B 


7.69 


8.25 


7.843E+08 


8.498E+07 


12091 


14699 


11294 


12997 


14250 


V 


7.69 


7.75 


2.068E+08 


5.365E+07 


12094 


14394 


11320 


12969 


13930 


M 


7.69 


7.25 


5.561E+07 


2.248E+07 


11973 


13558 


11281 


12676 


13070 


K 


7.69 


6.75 


1.542E+07 


6.641E+06 


11428 


12374 


10783 


11874 


11840 


Y 


7.39 


8.25 


8.636E+08 


8.379E+07 


12220 


15354 


11398 


13299 


14830 


AD 


7.39 


7.75 


2.191E+08 


5.455E+07 


12227 


14987 


11421 


13269 


14450 


P 


7.39 


7.25 


5.769E+08 


2.316E+07 


12120 


13985 


11392 


12952 


13430 


D 


7.39 


6.75 


1.575E+08 


6.754E+06 


11559 


12615 


10866 


12065 


12030 



photoionization models. However, in practice the different 
abundance techniques have different implicit assumptions 
about both the atomic physics relating to particular ionic 
species, and more generally, about the physics of the Hll 
region itself. Therefore, we should not expect the results to 
agree exactly. 

The physical and abundance parameters of the input 
models are given in Table [l] Note that we give three values 
of the ionization parameter, the nominal ionization param- 
eter input into the model, q, the ionization parameter at 
the inner surface of the ionized gas, qi n , and mean ioniza- 
tion parameter in the model, < ?(H) >. These are related 
to each other, but not in any obvious way. For example, q 
is calculated on the basis of an assumption that the tem- 



perature in the inner zone of the Hll region is 10,000 K. 
Because of our isobaric assumption, qi n is lower than q for 
the high-abundance models, and higher than qi n for the low- 
abundance models. The ratio of q ln to < q(R) > is always 
greater than two due to the spherical divergence of the radi- 
ation field, and can reach values as high as ten where radi- 
ation pressure gradients start to play a role in determining 
the gas pressure. 

As stated above, the input stellar spectrum, the physical 
conditions and the chemical abundance set were chosen so 
as to replicate, as far as p ossible, the observed sequence for 
extra galactic Hll regions (|Dopita et al.l 120001 ; iKewlev et al.l 
2006). In FigurcQJw e show the two main diagnostic diagrams 
from I Veilleux fc Osterbrockl II 19871 ) and the diagrams which 
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Figure 2. The relation between the oxygen abundance, 12 + log(0/H), and the electron temperature, T c given by the models in the 
high-ionization zone containing [Om] (left) and in the low- ionization zone containing the ions [Sn], [Oil], and [Nil] (right). The empirical 
fits we have derived for both high- and low-ionization zones as given by Eqs.[2]and[3]are shown as solid lines on these plots. 



iDopita et al.l (|2000f ) used to separate the chemical abun- 
dance from the ionization parameter at the high-abundance 
end of the scale. 

These models display quite a tight relationship between 
the oxygen abundance, 12+log(0/H), and the electron tem- 
perature, T a , given by the high-ionization zone from the 
[O ill] (A4959+A5007)/A4363 ratio, and in the low-ionization 
(H+ + He ) zone emitting the [On] (A3726+A3729), 
[Nil] A6584 and [Sn] (A6717+A6731) lines. This is shown 
in Fig [2] Note that there is larger scatter in the plot for the 
[O ill] zone, because the temperature in this zone is more 
sensitive to the input ionization parameter. On the left hand 
plot, the left-hand upper points correspond to high ioniza- 
tion parameter, while the reverse is true for the points in the 
lower right-hand side of the diagram. There is much less vari- 
ation in the local ionization parameter in the low-ionization 
zone of the model H n regions, so that the scatter is less for 
the right-hand panel of Fig. [2] 

For each zone of the model H n region we have generated 
empirical fits to the oxygen abundance - electron tempera- 
ture relationship. For [O ill] we find 

T e [Om] = 25200 - 1400a; - 2000 1 a: - 7.0| 17 , (2) 

while for the low-ionization zone we have 

T e (low) = 9500 + 400x - 2000|a; - 6.9| 235 ' (3) 

where x = 12 + log(0/H) in both cases. These curves em- 
phasize once again that the electron temperature is a key 
diagnostic for the chemical abundance in Hn regions. 



3 APPLICATION OF THE T c METHOD 

3.1 Physical conditions 

We have taken care to analyze the model spectra in exactly 
the same way as an observer would do when presented with 
the spectrum of real Hn regions. If there are no system- 
atic effects in the analysis technique, this exercise should 



yield the same electron temperatures and chemical abun- 
dances as was input into the models. One potential source 
of difference is the fact that the atomic data used in the 
nebular models is not necessarily the same as are used in 
the analysis of the spectra, although in many cases they 
are. The atomic data used for our models are detailed in 
ISutherland fc Dopital (1 19931 ) (see their Section 3.5.2), which 
are similar to the atomic data i ncluded in IRAF t ask te mden 
in the NEBULAR package by IShaw fc Dufouj l|l995h (see 
their Table 1). The output of the various nebular modeling 
codes have been compared against each other for a number 
of standard models on more than one occasion thr ough the 
so-ca lled Lexington Benchmarks code-off exercise (jFerlandl 
ll995T ). Any egregious errors made by any one of them have 
subsequently been eliminated from the appropriate codes. 
We can therefore be fairly c ertain that other ph otoionization 
codes such as Cloudy (e.g. lFerland et al.lll994T ) will produce 
very similar output to the Mappings Ills code used here. 

The model spectra were created without considering 
any extinction or reddening, and not stellar absorption in 
the H I Balmer lines, to decrease the uncertainties when an- 
alyzing the data. Thus, for the analysis of the model spectra 
we have used a reddening coefficient, c(H/3), and a equivalent 
width of the stellar absorption underlying the Hi Balmer 
lines, Wabs, equal to 0, following lLopez-Sanchez fc Estebanl 

(|20ogh - 

We derived the electron temperature, T e , and density, 
n e , of the ionized gas of the models using several emis- 
sion line ratios. The values obtained for each region are 
compiled in Table IA1I All determinations were computed 
with the IRAF task temden in the NEBULAR package 
l|Shaw fc Dufouiil9951 ). We used the updated atomic dataset 
for the 0+, S+, and S ++ ions as input to the NEBU- 
LAR routine. The refere nces are indicated in Table 4 of 
iGarcia-Roias et al.l (|2005l L 

In accord with Fig. [5] and with equations [2] and [31 we 
assumed a two-zone approximation for the Hn region to 
describe the temperature structure of the nebulae. In the 
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Figure 3. Comparison between the T c [Oni] given by the model (x-axis) and the temperature derived by analysis of the model spectra. 
The left hand panel directly compares the [O III] temperature from the model and the derived [O III] temperature, while the right -hand 
panel compares the fitted temperature drawn from the mean of the [O III] , [S III] and [Ar III] temperatures with the [O III] temperature 
from the model. The empirical fits we have derived are shown as a solid line for both panels. As in figure [2] the points are color coded 
according to the input abundance set. 



high-ionization H + + He + 
temperature, T e [Olll], and 
and T e [Arlll], T e (high), as 
for high ionization potential 
containing H + + He we 
T e [Sll] and T e [Nll], T e (lov/ 
tial ion temperatures. 

The line ratios used to 



zone, we used both the [O in] 
the mean of T e [Om], T e [Sm] 
the representative temperature 
ions. In the low-ionization zone 
adopted the mean of T e [O n] , 
), for the low ionization poten- 

measure the temperatures are: 



[Om] (A4959+A5007)/A4363, 

[S m] (A9069+A9530)/A6311, 

[Arm] (A7135+A7751)/A5192, 

[Nil] (A6548+A6583)/A5755, 

[On] (A3726+A3729)/(A7319+A7330), and 

[S n] (A6717+A6731)/(A4069+A4076). 

We note that some analyses of high-qualit y, deep spec- 
trophotometric data of bright H h gala xies (e.g. lHagele et al.l 
l2008l . l201ll ; |Perez-Montero et alj|20ich consider the electron 
temperature derived for a particular ion only to obtain the 
ionic abundance of that particular ion. However, the two- 
zone approximation for the Hll region is most widely used 
in extragalactic analyses. 

Since, in the analysis of the spectra, no information 
about the pressures or density adopted in the models was 
given, we used the standard density-sensitive line ratios 
[On] A3726/A3729, [Ni] A5198/A5200, [CI III] A5517/A5537 
and [S n] A6717/A6731 to derive n e . In practice, all models 
have a common pressure P/k = 3 x 10 5 cm -3 K, all density- 
sensitive ion line ratios should be at their low-density limit. 
Table |A"T1 compiles the value derived in all cases, and the av- 
erage value adopted for each model. Once n e was obtained 
we used it to derive T e using the [Om], [Sm], [Arm], [On], 
[S n], and [N n] line ratios, and we iterated until convergence 
was attained. 




4000 6000 8000 10000 12000 

T [low] (model) 

Figure 4. Comparison between the (emission weighted) T e for 
the low-ionization species [Nil], [On] and [Sn], as returned by 
the model (x-axis) and the T c derived for these same species from 
the temperature-sensitive line ratios (y-axis). The fit to the data 
(Eq. [Sj is also plotted. 



3.2 Results for T e 

The temperatures derived from the 

[Om] (A4959+A5007)/A4363 ratio, and the mean tem- 
perature derived for the high-ionization species [O ill] , [S ill] 
and [Ar m] , are compared with the [O ill] emission-line 
weighted temperature given by the model in Fig [3] The 
agreement for [O ill] is very good, while the mean of the 
high-ionization species tends to underestimate the true 
temperature. 

We empirically derive the following third-order correc- 
tion formula to transform the measured [Om] line temper- 
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Table 2. Total abundances derived from the analysis of line intensity ratios of the model spectra using IRAF task ionic considering 
the electron temperature derived from all the high-ionization species, T (high), and the low ionization species, T G (low), separately. All 
abundances are expressed in the form of 12 + log (X/H) or log (X/O). 



Model 


O++/O+ 


O/H 


N/H 


N/O 


S/H 


S/O 


Ne/H 


Ne/O 


Ar/H 


Ar/O 


Cl/H 


Cl/O 


A 


-0.592 


8.816 


7.432 


-1.384 


7.147 


-1.669 


8.374 


-0.443 


6.883 


-1.933 


4.424 


-4.393 


B 


0.326 


7.777 


6.129 


-1.648 


5.974 


-1.803 


7.109 


-0.668 


5.512 


-2.264 


3.426 


-4.351 


C" 


-1.771 


8.880 


8.264 


-0.616 


7.601 


-1.279 


9.125 


0.245 






4.814 


-4.066 


D 


-0.680 


7.571 


5.923 


-1.649 


5.843 


-1.728 


7.230 


-0.342 


5.618 


-1.954 


3.240 


-4.332 


E 


-1.258 


8.197 


6.629 


-1.567 


6.466 


-1.731 


8.165 


-0.032 


6.632 


-1.565 


3.761 


-4.436 


F 


-0.126 


8.495 


6.959 


-1.536 


6.810 


-1.686 


7.912 


-0.583 


6.340 


-2.155 


4.124 


-4.371 


G 


-0.774 


9.157 


8.105 


-1.051 


7.595 


-1.562 


8.644 


-0.512 


7.301 


-1.855 


4.822 


-4.335 


II" 


-0.889 


8.763 


8.205 


-0.558 


7.669 


-1.094 


8.600 


-0.163 


7.441 


-1.323 


4.884 


-3.879 


I 


-0.654 


9.050 


7.796 


-1.254 


7.306 


-1.744 


8.495 


-0.555 


7.033 


-2.017 


4.568 


-4.482 


J" 


-0.267 


8.560 


8.080 


-0.480 


7.723 


-0.837 


8.164 


-0.396 


7.011 


-1.548 


4.905 


-3.654 


K 


-0.697 


7.857 


6.235 


-1.623 


6.142 


-1.715 


7.530 


-0.327 


5.931 


-1.927 


3.444 


-4.413 


L 


-0.925 


8.362 


6.807 


-1.555 


6.627 


-1.735 


8.103 


-0.259 


6.573 


-1.789 


3.914 


-4.448 


M 


-0.067 


7.805 


6.153 


-1.653 


6.094 


-1.711 


7.212 


-0.594 


5.597 


-2.208 


3.423 


-4.383 


N 


-0.103 


8.733 


7.332 


-1.401 


7.116 


-1.617 


8.160 


-0.573 


6.630 


-2.103 


4.411 


-4.323 


O 


-1.863 


9.358 


8.320 


-1.038 


7.611 


-1.747 


9.332 


-0.026 


7.996 


-1.362 


4.835 


-4.523 


P 


-0.036 


7.502 


5.831 


-1.670 


5.786 


-1.716 


6.900 


-0.602 


5.283 


-2.219 


3.121 


-4.381 


Q 


0.236 


8.079 


6.456 


-1.622 


6.308 


-1.770 


7.422 


-0.656 


5.829 


-2.250 


3.721 


-4.358 


R 


-0.647 


9.073 


7.877 


-1.196 


7.439 


-1.633 


8.597 


-0.476 


7.174 


-1.899 


4.698 


-4.374 


S 


-0.255 


8.302 


6.708 


-1.594 


6.598 


-1.704 


7.747 


-0.555 


6.163 


-2.139 


3.900 


-4.402 


T 


0.112 


8.256 


6.657 


-1.599 


6.539 


-1.717 


7.621 


-0.635 


6.036 


-2.221 


3.896 


-4.360 


U 


0.160 


8.456 


6.920 


-1.536 


6.760 


-1.696 


7.824 


-0.632 


6.252 


-2.204 


4.117 


-4.338 


V 


0.230 


7.780 


6.129 


-1.651 


6.010 


-1.770 


7.125 


-0.655 


5.524 


-2.256 


3.423 


-4.357 


w 


-0.677 


8.564 


7.055 


-1.510 


6.847 


-1.717 


8.167 


-0.397 


6.634 


-1.930 


4.132 


-4.432 


x° 


-0.567 


8.303 


7.803 


-0.501 


7.492 


-0.811 


7.765 


-0.538 


6.544 


-1.759 


4.640 


-3.663 


Y 


0.309 


7.486 


5.819 


-1.667 


5.664 


-1.822 


6.814 


-0.672 


5.211 


-2.275 


3.127 


-4.359 


z 


-2.018 


9.183 


8.028 


-1.155 


7.455 


-1.728 


9.371 


0.188 


7.939 


-1.244 


4.695 


-4.488 


AA 


-1.189 


9.320 


8.267 


-1.053 


7.622 


-1.698 


8.962 


-0.357 


7.644 


-1.676 


4.854 


-4.466 


AB 


-0.447 


8.149 


6.537 


-1.612 


6.429 


-1.720 


7.672 


-0.477 


6.077 


-2.072 


3.728 


-4.422 


AC 


-1.232 


9.163 


7.988 


-1.174 


7.449 


-1.714 


8.923 


-0.239 


7.563 


-1.600 


4.699 


-4.464 


AD 


0.243 


7.481 


5.813 


-1.668 


5.697 


-1.784 


6.822 


-0.660 


5.218 


-2.264 


3.123 


-4.358 


AE 


0.048 


8.092 


6.467 


-1.625 


6.373 


-1.719 


7.466 


-0.625 


5.868 


-2.224 


3.719 


-4.373 


AF 


-1.473 


8.599 


7.132 


-1.467 


6.870 


-1.729 


8.644 


0.045 


7.177 


-1.423 


4.151 


-4.448 


AG 


-1.281 


8.873 


7.525 


-1.348 


7.161 


-1.713 


8.759 


-0.114 


7.347 


-1.526 


4.429 


-4.444 



a For these models, there is not estimation of the high and low electron temperatures using the [O III] and [Oil] ions because their auroral lines are not 
tabulated (sec Table fA It . Furthermore, ionic task does not allow electron temperatures lower than 5,000K for the low ionization ions, so we adopted 
T c (low) = 5,000 K in these four cases. 



attires into the model [Oiii] line temperatures: 
it =t 4 (l-0.05(i 4 - l) 2 ), 



(4) 



where ti = T e /10,000K is the temperature measured from 
the spectra, and t% is the temperature corrected to that de- 
livered by the corresponding model in units of 10 4 K. The 
corrections to T c [Oiii] are almost always smaller than the 
intrinsic uncertainty in the derived temperature - typically 
of order 500-1000 K. 

For the mean of the high excitation species [O ill] , [S ill] 
and [Arm] the corresponding correction formula is: 



t\ =£4(1.12 



0.02*1) 



0.087, 



(•») 



with *4 as the average value of the electron temperature 
derived using the [O m] , [S m] and [Ar ill] ratio. 

For the low- ionization zones of the model Hll regions 
we use the mean of the [0 11] , [S 11] and [N 11] temperatures 
given by the models, as was also done in the fitting exercise. 
The result of the comparison is shown in Fig [3] As we see, 
we find a clear linear correlation between model and fit, 
although the slope is not exactly unity. The empirical fit is 
given by the formula: 



t% = t 4 (l-28 - 0.03*4) - 0.15, 



(6) 



where £4 and t% are the temperature measured from the spec- 
trum and the corrected to that delivered by the correspond- 
ing model in units of 10 4 K, respectively. The correlation 



coefficients of the fits in Eqn.|31[S]and[S]are all in excess of 
0.98. 

These differences between the temperatures delivered 
by the model are the result of the temperature struc- 
ture within the ionized region. The temperatures delivered 
by the Mappings Ills code are weighted by the emissiv- 
ity of, for example, the [Om] A5007 line for [Om]. The 
[O m] A4363 line is more stron gly affected by temperatur e 
stratification as pointed out bv lPeimbert fe Costerd IJ1969), 
so that the temperature delivered by measurement of the 
[Om] (A4959+A5007)/A4363 ratio is different to the strong- 
line weighted mean temperature delivered by the code. 

For observers who are only able to measure the [O m] 
temperature in the high-ionization zone (H + + He + ), and 
need a correction formula to estimate the temperature in 
the low-ionization (H + + He ) zone, the models provide a 
convenient relationship that fits to an accuracy of ±300 K 
for temperatures above 10,000 K: 

Te(low) = T e [0 m]+450-70exp[(T e [O m]/5, 000) 1 ' 22 ]. (7) 

Note that the temperature sensitivity to the abundance 
is strong, and that in high-abundance Hll regions, the 
low ionization species deliver a lower electron tempera- 
ture. Note however, that the form and s hape of Eq. [7] i s 
rather different from the one proposed by iGarnetj l|1992h . 
T e (low) = 0.7 x T e [0 ill] + 3, 000. However, both the tem- 
perature offset and the slope are similar in the range 
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Figure 5. A comparison of the total abundances derived from the analysis of line intensity ratios of the model spectra and using T c (high) 
and T c (low) as representative for the electron temperature of the high- and low-ionization species, respectively (see Tabled, with the 
total abundances provided by the models. The fit is fairly good for the species for which the optical line emission arises principally in the 
low-ionization zone, N, S and CI, but there is clearly a systematic error for the other species, O, Ne and Ar. For 12+log(0/H) = 9.39, 
the T e method fails because IRAF task ionic does not allow electron temperatures lower than 5,000 K for the low ionization ions, and 
hence T c (low) = 5,000 K was adopted. 
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Figure 6. A comparison of the total abundances derived by using T e [Olll] using the [Om] (A4959+A5007)/A4363 ratio in the high- 
ionization zone and the T e [Oll] using the [On] (A3726+A3729)/(A7319+A7330) ratio in the low-ionization zone (see Table [3]l with the 
total abundances provided by the models. As in Fig. \E\ the fit is fairly good for the species for which the optical line emission arises 
principally in the low-ionization zone, N, S and CI. Although it is slightly overprcdicting the model abundances, the fit is now acceptable 
for O, except for 12+log(0/H) = 9.39 because of the aforementioned problem when calculating T e [Oll] for the low ionization ions. 
However, the fit still has a systematic error in the case of Ne (which is seen only in its high-ionization stage Ne ++ and Ar (for which the 
emission is also dominated by its high ionization lines, and which is strongly affected by charge-exchange reactions). Random errors are 
also larger for these species. 
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Table 3. Total abundances derived from the analysis of line intensity ratios of the model spectra using IRAF task ionic considering 
only the electron temperature derived from the [Oin] lines in the high-ionization zone, T e [Oin], and that derived from the [On] lines, 
T c [On], in the low-ionization zone. All abundances are expressed in the form of 12 + log (X/H) or log (X/O). 



Model 


O++/O+ 


O/H 


N/H 


N/O 


S/H 


S/O 


Ne/H 


Ne/O 


Ar/H 


Ar/O 


Cl/H 


Cl/O 


A 


-0.511 


8.707 


7.362 


-1.344 


7.114 


-1.593 


8.269 


-0.438 


6.797 


-1.910 


4.388 


-4.319 


B 


0.344 


7.616 


6.021 


-1.595 


5.884 


-1.732 


6.968 


-0.648 


5.406 


-2.210 


3.337 


-4.279 


C" 


-1.771 


8.880 


8.264 


-0.616 


7.601 


-1.279 


9.125 


0.245 






4.814 


-4.066 


D 


-0.662 


7.435 


5.828 


-1.607 


5.770 


-1.666 


7.109 


-0.327 


5.521 


-1.915 


3.071 


-4.365 


E 


-1.239 


8.065 


6.538 


-1.527 


6.397 


-1.669 


8.047 


-0.018 


6.538 


-1.527 


3.686 


-4.380 


F 


-0.0S7 


8.397 


6.896 


-1.501 


6.759 


-1.638 


7.823 


-0.574 


6.271 


-2.126 


4.076 


-4.321 


G 


-0.559 


9.027 


8.031 


-0.996 


7.609 


-1.418 


8.511 


-0.515 


7.185 


-1.842 


4.837 


-4.190 


II' 1 


-0.934 


8.758 


8.200 


-0.558 


7.649 


-1.109 


8.598 


-0.160 


7.443 


-1.315 


4.860 


-3.898 


I 


-0.140 


8.888 


7.730 


-1.158 


7.408 


-1.480 


8.322 


-0.566 


6.891 


-1.997 


4.671 


-4.217 


J a 


-0.267 


8.560 


8.080 


-0.480 


7.723 


-0.837 


8.164 


-0.396 


7.011 


-1.548 


4.905 


-3.654 


K 


-0.698 


7.740 


6.153 


-1.587 


6.073 


-1.667 


7.427 


-0.313 


5.849 


-1.891 


3.371 


-4.369 


L 


-0.903 


8.236 


6.720 


-1.516 


6.562 


-1.674 


7.990 


-0.247 


6.482 


-1.755 


3.844 


-4.392 


M 


-0.043 


7.675 


6.064 


-1.610 


6.024 


-1.651 


7.095 


-0.579 


5.506 


-2.168 


3.351 


-4.323 


N 


-0.010 


8.654 


7.291 


-1.363 


7.100 


-1.554 


8.084 


-0.570 


6.575 


-2.079 


4.388 


-4.266 


O 


-1.759 


9.192 


8.208 


-0.984 


7.549 


-1.643 


9.172 


-0.020 


7.887 


-1.305 


4.768 


-4.424 


P 


-0.013 


7.367 


5.740 


-1.627 


5.713 


-1.654 


6.780 


-0.586 


5.189 


-2.177 


3.047 


-4.320 


Q 


0.248 


7.943 


6.364 


-1.579 


6.232 


-1.711 


7.303 


-0.640 


5.737 


-2.206 


3.645 


-4.298 


R 


-0.508 


8.962 


7.811 


-1.151 


7.429 


-1.533 


8.486 


-0.476 


7.080 


-1.882 


4.686 


-4.276 


S 


-0.228 


8.185 


6.629 


-1.556 


6.538 


-1.647 


7.642 


-0.543 


6.080 


-2.105 


3.837 


-4.348 


T 


0.139 


8.143 


6.583 


-1.560 


6.480 


-1.663 


7.520 


-0.623 


5.958 


-2.185 


3.835 


-4.308 


U 


0.198 


8.363 


6.862 


-1.501 


6.717 


-1.646 


7.740 


-0.623 


6.189 


-2.174 


4.071 


-4.292 


V 


0.243 


7.639 


6.033 


-1.606 


5.930 


-1.709 


7.001 


-0.638 


5.428 


-2.210 


3.344 


-4.294 


W 


-0.645 


8.451 


6.977 


-1.473 


6.793 


-1.658 


8.063 


-0.388 


6.550 


-1.901 


4.074 


-4.377 


X° 


-0.567 


8.303 


7.803 


-0.501 


7.492 


-0.811 


7.765 


-0.538 


6.544 


-1.759 


4.640 


-3.663 


Y 


0.297 


7.319 


5.701 


-1.619 


5.562 


-1.757 


6.670 


-0.650 


5.098 


-2.221 


3.031 


-4.289 


Z 


-1.929 


9.027 


7.922 


-1.105 


7.393 


-1.634 


9.222 


0.195 


7.844 


-1.183 


4.626 


-4.400 


AA 


-1.054 


9.173 


8.172 


-1.002 


7.591 


-1.583 


8.818 


-0.356 


7.522 


-1.651 


4.821 


-4.352 


AB 


-0.425 


8.022 


6.450 


-1.572 


6.362 


-1.660 


7.558 


-0.464 


5.986 


-2.036 


3.657 


-4.365 


AC 


-1.134 


9.022 


7.895 


-1.127 


7.406 


-1.616 


8.787 


-0.235 


7.450 


-1.571 


4.652 


-4.370 


AD 


0.251 


7.330 


5.709 


-1.621 


5.610 


-1.720 


6.690 


-0.641 


5.116 


-2.215 


3.037 


-4.293 


AE 


0.071 


7.967 


6.383 


-1.584 


6.306 


-1.661 


7.355 


-0.611 


5.781 


-2.185 


3.650 


-4.316 


AF 


-1.445 


8.476 


7.046 


-1.429 


6.809 


-1.667 


8.531 


0.056 


7.088 


-1.388 


4.083 


-4.392 


AG 


-1.214 


8.740 


7.434 


-1.305 


7.109 


-1.630 


8.633 


-0.106 


7.244 


-1.495 


4.373 


-4.367 



a For these models there is not estimation of the high and low electron temperatures using the [O III] and [Oil] ions because their auroral lines arc not 
tabulated (sec Table fA It . The results shown here are the same that in Tablc[2]bccausc we used the T c (high) and T c (low). Furthermore, ionic task does not 
allow electron temperatures lower than 5,000 K for the low ionization ions, so we adopted T c (low) — 5,000 K in these four cases. 



8,000 ^ T e [0 in] sC 12,000 where most of the high-quality 
nebula spectra have been obtained. 



3.3 Ionic abundances 

The IRAF package NEBULAR (|Shaw fe Dufourlll995l ) has 
been used to analyze the intensities of the forbidden lines 
using the temperatures and densities derived for the high- 
and low-ionization zones to derive ionic abundances of + , 
0++, N+, S+, S++, Ne++, Ar++, Ar+ 3 , C1++, and C1+ 3 . 
The ionic abundances so derived are listed in the Appendix 
in Tabled 



3.4 Total abundances 

In all cases we adopted O/H = 0+/H+ + 0++/H+ to de- 
termine the total oxygen abundance. Although a very weak 
nebular Hen A4686 line is produced in several models, the 
relative contribution of He ++ to the total amount of helium 
is negligible, implying that 3+ has also a very low abun- 
dance in the nebula, thus we did not consider its contribution 
to the total O/H ratio. For the rest of the elements, we have 
to adopt a set of ionization correction factors (ICFs) to cor- 
rect for the unseen ionization sta ges. The ICFs adopted here 
are b asically the same used by lLopez-Sanchez fc Estebanl 
J2009I ). 



To deriv e the nitrogen abundance we assumed the stan- 
dard ICF bv lPeimbert fc Costerd (|l969T l: N/O = N+/0+, 
which is the typical assumption considered in the anal- 
ysis of the ionized gas in extr agalactic Hn regions (e.g. 
lLopez-Sanchez fc Estebanll20ld '). 

We measured two ionization stages of sulphur, S + and 
S ++ , in all spectra. However, a significant contr ibution of 
S 3+ is expected. We adopted the ICF given by IStasiriskal 
i|1978T) , which is based on photoionization models of H II re- 
gions and is expressed as a function of the + /0 ratio. 

In the case of neon we also applie d the classical ICF 
proposed bv lPeimbert fc Costerd (1 19691 ), that assumes that 
the ionization structure of Ne is similar to that of O. This 
is a good approximation for high ionization objects, where 
a small fraction of Ne + is expected. 

For argon we have determinations of the Ar ++ and Ar 3+ 
abundances. However, some contribution of Ar + is expected. 
Hence, the total argon abun dance was ca l culate d by consid- 
ering the ICFs proposed bv llzotov et all (| 19941 ). 

We measured lines Cl ++ and Cl 3+ . As can be seen in 
Table IA2I the dominant ionization stage is Cl ++ and the 
contribution of Cl + to the total abundance is rather small. 
To take int o account the Cl + fraction we have adopted the 
relation bv lPeimbert fc Torres-Peimbertj (|l977T ). 

The results we obtain depend on whether we use 
all of the high-ionization species to determine T e in the 
high-ionization zone, T e (high), and all the low-ionization 
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species to compute T e in the low- ionization zone, T e (low), 
or whether we use only the temperature derived from [O ill] 
line ratio, T c [O ill] . The total abundances we obtain by these 
two approaches are shown in Tables [2] and [3] and are shown 
in graphical form in Fig [5] and [6] respectively. 

When all of the high ionization species are used 
to derive T c (high), the oxygen abundance is systemati- 
cally overestimated by ~ 0.2 dex and the other species 
seen principally in their high-excitation stages, Ne and 
Ar, are overestimated by somewhat larger amounts. For 
12+log(0/H)=9.39, the T c method fails because IRAF 
task ionic does not allow electron temperatures lower than 
5,000 K for the low ionization ions. If we use only the temper- 
ature derived from the [O ill] lines, T e [O ill] using the [O ill] 
(A4959+A5007)/A4363 ratio in the high-ionization zone, and 
the temperature derived from the [On] lines, T e [On] using 
the [On] (A3726+A3729)/(A7319+A7330) ratio in the low- 
ionization zone, the systematic error is much less, specially 
for the total oxygen abundance, although it is still slightly 
overpredicted. 

Although there are small systematic errors in this 
procedure, all are correctable. It is clear that the T e 
method returns reliable abundances. Provided that the 
[Om] (A4959+A5007)/A4363 ratio is used to determine the 
T e in the high-ionization zone, the scatter in the derived 
abundances of O, N, S and CI is typically ±0.1 dex, provided 
that the O abundance is below 12+log(0/H) < 9. For high 
abundances, the T a method returns systematically low abun- 
dances, wh ere it can be app lied. This is in agreement with 
the work of lStasiriskal IJ2002r) and lStasih ska (2005), who pre- 
dicted that temperature gradients in these high-metallicity 
H II regions can cause the abundances to be underestimated 
by as much as ~ 0.4 dex. 

Both for Ne and Ar, the systemic error and the scatter 
are distressingly large. We do not think that this is a con- 
sequence of using slightly different atomic data for models 
and line analysis using IRAF. However the cause is likely 
to be different in the two cases. In the case of Ne, the error 
is probably caused by the fact that only Ne ++ is observ- 
able, and the [Nelll] AA3869, 3968 doublet has sufficiently 
high excitation energy that the line emissivity is strongly 
biased towards the hottest part of the Hll region. In the 
case of Ar, charge-exchange processes are very important in 
determining the ionic balance, and these strongly affect the 
ionizati on balance and the refore the ionization correction 
factors (jDopita et al.lll997l ). Furthermore, the assumptions 
of the ICFs used to derived the total abundances of Ne and 
Ar from t heir ionic abundances may be not valid in some 
cases (e.g. iPerez-Montero et aUl2007n . 

Given that Fig [2] shows that the electron temperature 
and the abundance are closely correlated, it is interesting 
to see just how well the oxygen abundance can be recov- 
ered from the measured [Om] temperature alone. To do 
this we used our empirical relation (Eq. Q to correct the 
measured temperature and then applied the relationship be- 
tween T e [Om] and 12+log(0/H) given by Eq. [2] to obtain 
the oxygen abundance. The result is shown in Fig. [7] For 
abundances 12+log(0/H) > 8.0 the abundances are recov- 
ered to an accuracy of ±0.1 dex. At lower abundances, sig- 
nificant scatter is introduced by variations in the ionization 
parameter q, and it becomes essential to determine this pa- 
rameter independently. 




8.0 8.5 9.0 

12 + log(0/H) (model) 

Figure 7. A comparison of the abundances derived by using only 
the T e [Oni] measured from the [Om] (A4959+A5007)/A4363 ra- 
tio. For abundances 12+log(0/H) > 8.0 the abundances are re- 
covered to an accuracy of ±0.1 dex. At lower abundances, the 
scatter is much increased because of the increased sensitivity of 
T e [Oni] with the ionization parameter in this region. 



4 APPLICATION OF THE SEL METHOD 

The strong emission line (SEL) technique to derive nebular 
abundances is typically used when the spectra are of too 
low signal-to-noise ratio for the temperature sensitive emis- 
sion lines to be detectable. This typically encompasses the 
majority of global determinations of metallicity in external 
galaxies, and all galaxies with redshifts large enough to be 
of cosmological interest. An extensive review of 10 metallic- 
ity calibrations, including th eoretical and emp i rical m eth- 
ods, has been presented by iKewle v & Ellison (200 8j), us- 
ing data gathered from the SDSS. lLopez-Sanchez fc Estebanl 
(2010) also review the most common empirical calibrations 
and compare their results with those derived using the T 
method. The majority of the empirical calibrations rely on 
ratios between bright emission lines to estimate the oxygen 
abundance. The most commonly used of these are known by 
their shorthand contractions; R23, S23, P, N2, O3N2, N2O2 
and y. Their definitions are: 



R» 



R2 = 



I([Q m])A4959 + J([Q iii])A5007 
H/3 

J([On])AA3726,3729 



H/3 



R23 =R 3 + R2, 
R 3 



P = 



R23' 



, #3 . 1 

y = log _ = ]og __ 



N 2 = 



J([Nn])A6584 
Ha 



(8) 

(9) 

(10) 

(11) 
(12) 

(13) 
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Table 4. Results for the oxygen abundance, in the form 12+log(0/H), according to the most-commonly used empirical SEL calibrations. 



Model 


Tc 


BRANCH 


M91 


KD02o 


KK04 


KN202 


P01 


PT05 


PVT10 


PM11 


D02 


PP04a 


PP04b 


PP04c 


S23 


S23023 


A 


8.71 


HIGH 


8.53 


8.84 


8.67 


8.64 


8.21 


8.11 


8.44 


8.49 


8.56 


8.46 


8.41 


8.40 


8.96 


8.60 


B 


7.62 


LOW 


7.50 


7.67 


7.81 




7.39 


7.36 


7.17 


7.12 


7.62 


7.73 


7.74 


7.94 


7.45 


7.74 


C 


8.88 


HIGH 


9.13 


9.31 


9.35 


9.24 


8.86 


8.43 


9.11 


8.94 


8.72 


8.59 


8.58 


9.25 


8.53 


9.65 


D 


7.44 


LOW 


7.23 


7.57 


7.53 




7.34 


7.01 


6.73 


6.31 


7.71 


7.80 


7.84 


8.28 


7.28 


8.31 


E 


8.01 


LOW 


7.96 


7.85 


8.10 




8.77 


6.94 


7.18 


6.74 


8.21 


8.19 


8.18 


8.51 


8.17 


8.50 


Fl" 


8.40 


LOW 


8.17 


8.28 


8.36 




8.02 


8.13 


7.93 


7.95 


8.26 


8.23 


8.20 


8.17 


8.74 


8.37 


Fu" 


8.40 


HIGH 


8.17 


8.73 


8.55 


8.34 


8.27 


8.22 


7.93 


7.95 


8.26 


8.23 


8.20 


8.17 


8.74 


8.37 


G 


9.03 


HIGH 


9.07 


9.22 


9.16 


9.09 


8.84 


8.70 


8.73 


8.76 


8.54 


8.45 


8.39 


8.69 


8.46 


9.30 


II 


8.76 


HIGH 


9.14 


9.33 


9.35 


9.26 


8.88 


8.57 


8.94 


8.88 


8.66 


8.54 


8.51 


9.05 


8.39 


9.67 


I 


8.89 


HIGH 


8.95 


9.02 


9.06 


8.92 


8.76 


8.69 


8.65 


8.67 


8.44 


8.37 


8.31 


8.45 


8.55 


8.94 


J 


8.56 


HIGH 


9.22 


9.37 


9.47 


9.28 


8.91 


8.69 


9.01 


9.00 


8.49 


8.41 


8.35 


8.98 


8.04 


9.74 


K 


7.74 


LOW 


7.54 


7.81 


7.79 




7.78 


7.15 


7.15 


6.74 


7.94 


7.98 


8.03 


8.31 


7.67 


8.33 


L 


8.24 


LOW 


8.12 


7.97 


8.27 




8.57 


7.52 


7.67 


7.28 


8.32 


8.28 


8.24 


8,41 


8.42 


8.42 


M 


7.68 


LOW 


7.50 


7.58 


7.80 




7.40 


7.41 


7.18 


7.07 


7.78 


7.85 


7.91 


8.07 


7.62 


8.07 


N 


8.65 


HIGH 


8.61 


8.77 


8.74 


8.63 


8.41 


8.36 


8.52 


8.55 


8.38 


8.32 


8.27 


8.26 


8.74 


8.51 


O 


9.19 


HIGH 


9.00 


9.07 


9.14 


9.06 


8.65 


8.29 


8.76 


8.67 


8.81 


8.65 


8.70 


8.97 


9.02 


9.29 


P 


7.37 


LOW 


7.1(1 


7.53 


7.52 




7.00 


7.20 


6.75 


6.62 


7.54 


7.67 


7.64 


8.05 


7.24 


8.04 


Q 


7.94 


LOW 


7.01 


7.98 


8.15 




7.73 


7.74 


7.53 


7.52 


7.86 


7.92 


7.98 


7.98 


7.90 


7.86 


R 


8.96 


HIGH 


8.88 


8.97 


9.02 


8.92 


8.61 


8.51 


8.56 


8.57 


8.61 


8.51 


8.46 


8.54 


8.86 


8.96 


S 


8.19 


LOW 


8.07 


8.13 


8.27 




8.00 


7.98 


7.77 


7.64 


8.17 


8.16 


8.16 


8.17 


8.34 


8.22 


T 


8.14 


LOW 


8.05 


8.08 


8.26 




7.85 


7.92 


7.67 


7.68 


8.02 


8.04 


8.08 


8.05 


8.20 


8.04 


Ul° 


8.36 


LOW 


8.13 


8.22 


8.32 




7.90 


8.02 


7.88 


7.89 


8.13 


8.12 


8.13 


8.08 


8.40 


8.13 


Un° 


8.36 


HIGH 


8.13 


8.58 


8.54 


8.33 


8.34 


8.29 


7.88 


7.89 


8.13 


8.12 


8.13 


8.08 


8.40 


8.13 


V 


7.64 


LOW 


7.50 


7.65 


7.81 




7.38 


7.36 


7.16 


7.11 


7.65 


7.75 


7.78 


7.97 


7.51 


7.83 


Wl a 


8.45 


LOW 


8.27 


8.29 


8.41 




8.45 


8.04 


7.99 


8.00 


8.43 


8.36 


8.31 


8.35 


8.72 


8.42 


Wu a 


8.45 


HIGH 


8.27 


8.78 


8.59 


8.38 


8.11 


7.99 


7.99 


8.00 


8.43 


8.36 


8.31 


8.35 


8.72 


8.42 


X 


8.30 


HIGH 


9.32 


9.40 




9.28 


8.92 


8.62 


9.12 


9.18 


8.36 


8.31 


8.26 


9.05 


7.71 


9.76 


Y 


7.32 


LOW 


7.15 


7.47 


7.52 




7.01 


7.13 


6.76 


6.68 


7.39 


7.55 


7.38 


7.92 


7.12 


7.70 


z 


9.03 


HIGH 


8.76 


9.01 


8.98 


8.92 


8.36 


7.98 


8.73 


8.61 


8.82 


8.66 


8.72 


8.95 


9.29 


9.08 


AA 


9.17 


HIGH 


9.03 


9.17 


9.14 


9.07 


8.73 


8.52 


8.66 


8.65 


8.71 


8.58 


8.57 


8.77 


8.81 


9.28 


AB 


8.02 


LOW 


7.90 


7.84 


8.12 




7.96 


7.09 


7.61 


7.36 


8.11 


8.11 


8.13 


8.22 


8.12 


8.25 


AC 


9.02 


HIGH 


8.79 


9.02 


8.97 


8.92 


8.42 


8.18 


8.53 


8.53 


8.77 


8.63 


8.65 


8.70 


9.12 


9.00 


AD 


7.33 


LOW 


7.14 


7.46 


7.52 




6.99 


7.14 


6.75 


6.66 


7.42 


7.57 


7.43 


7.94 


7.15 


7.79 


AE 


7.97 


LOW 


7.90 


7.89 


8.14 




7.75 


7.70 


7.55 


7.51 


7.94 


7.98 


8.03 


8.05 


8.02 


8.02 


AFP 


8.48 


LOW 


8.38 


7.85 


8.44 




9.43 


7.29 


7.80 


7.93 


8.51 


8.42 


8.37 


8.62 


8.84 


8.61 


AFu a 


8.48 


HIGH 


8.38 


8.80 


8.72 


8.47 


8.07 


7.75 


7.80 


7.93 


8.51 


8.42 


8.37 


8.62 


8.84 


8.61 


AG 


8.74 


HIGH 


8.45 


8.76 


8.65 


8.67 


8.03 


7.74 


8.48 


8.50 


8.69 


8.56 


8.54 


8.62 


9.15 


8.71 



a For these models, for which wc derived 12+log(0/H)'^8.4 following the T c method, we list the results of the empirical calibrations considering both the 

low and the high mctallicity branches. 

NOTE: The empirical calibrations and the parameters used for each of them arc: 

M91: iMcGaugh ( 199J1) using J? 2 3 and y; 

KD02: Kcwlcy &c Dopita (2002) using P23 and the ionization parameter defined in that paper, <lKDQ2o'i 

KK04: Kobulnicky &c Kcwlcy (2004) using P23 and the ionization parameter defined in that paper, qKKO&'i 

KN202: Kcwlcy &c Dopita (2002) using the N 2 2 parameter (calibration only valid for objects in the high mctallicity branch); 

P01: Pilyugin (2001a,b) using R 23 and P: 

PT05: [Pilvugin fc T huan (200(1) using K23 and P: 

PVT10: Pilyugin ct al. (2010) using R 23 , P, N 2 and S 2 ; 

PM11: Pilyugin & Mattsson I 20"Tl|) using R 3 , N 2 and S 2 ; 

D02: Denicolo, Terlcvich & Tcrlevich (2002) using the N 2 parameter; 

PP04a: Pcttini Sz Pagcl (2004), using a linear fit to the N 2 parameter; 

PP04b: Pettini &; Pagel (2004), using a cubic fit to the N 2 parameter; 

PP0 4c: Pettini &: Pagel (2004), us ing the 3 N 2 parameter; 

S23: Pcrcz-Montcro & Diaz (2005) using S 23 ; 

S23023: Pcrcz-Mont cro & Di'aj J2005) using S 2 3 and R 2 3- 



O3N2 



N2O2 



S23 



_ I([Oin])A5007 

~ I([Nn])A6584' 

J([Nn])A6584 
~ J([Oii])AA3726,3729' 

J([S n])AA6717, 6731 + J([S m])AA9076, 9532 



H/3 



023(^23 = ■= — • 
-K23 



Si mple line ratio s 



(14) 
(15) 
(16) 
(17) 



such as J?23 (I Jensen et all 
&3 fVi'lchcz fc Estebanl 1 19961 ; 



19761: iPaeel et all I1979T ) 

Christensen et all Il997l ; iDi'az fc Perez-Monterol 120001) 

indeed any forbidden line to recombination line ratio- suffer 



or 



from being two- valued, and hence a separate calibration 
for the low and high metallicity regimes is usually needed. 
In the case of the R23 index, the calibrations are given for 
12+log(0/H)<8.1 (low metallicity) and 12+log(0/H)>8.4 
(high metallicity). That means that a very large fraction 
of the star-forming regions lie in the ill-defined turning 
zone around 12+log(O/H)~8.20-8.30, where regions with 
the same -R23 value have oxygen abundances that differ by 
al most an order of magnitude (see, f or example, Figure A.l 
in lLopez-Sanchez fc Estebanl 120101 ). The reason of this 
behavior is that the intensity of the oxygen (or any other 
heavy element) lines do not monotonically increase with 
metallicity. At low abundance, the lines of the heavy 
element being considered are weak relative to a hydrogen 
recombination line due to their low abundance. However, 
at high abundance they are weak due to the low electron 



14 A.R. Lopez-Sanchez et al. 



temperature, which suppresses collisional excitations into 
the excited state responsible for the forbidden line. The 
turn-around in the line ratio can be extended to higher 
abundance by choosing lines with lower excitation poten- 
tials which make the ratio less temperature-sensitive. Hence 
the utility of 5*23 over 7?23- The use of multiple line ratios 

also helps to remove the abundance ambiguities. 

Ratios such as C^A ^Alloin et al.|[l979l ; |Pettini fc Pagell 

|200J) were introduced in an attempt to sidestep the abun- 
dance ambiguity entirely. This ratio is much more mono- 
tonic in abundance because of the strong secondary compo- 
nent of N enrichment alluded to above; see Eq.Q] This helps 
at high abundance, but the ratio once again becomes am- 
biguous when 12 + log(Q/H) < 8.0. A rati o such as N2O2 
JDopita et alj|200d : IKewlev fc Dopitall2002l ) also provides a 
useful calibration for metallicity at high abundance. Its ad- 
vantage is that both N + and + co-exist in the same zone 
of the nebula, the [O ill] lines become weaker at high metal- 
licity due to the low electron temperature, while the [Nil] 
lines become stronger due to the higher relative abundance 
of this element. However, its disadvantage is that the ni- 
trogen and oxygen emission lines are widely separated in 
wavelength, making flux calibration and reddening correc- 
tions more critical. 

The fundamental weakness in the use of all these strong- 
line ratios is that the overall spectrum of an H 11 region 
depends not only upon the chemical abundances, but also 
upon the ionization parameter q or U. That this is so is 
abundantly o bvious from Fig . 1. Th e importance of this was 
recognized bv lBaldwin et al.l |l98lj), who were amongst the 
first to use line ratio diagnostic diagrams. However, the first 
empirical cal ibration involvin g an ionization parameter was 
presented bv lMcGaugbl (|l99ll ). who developed mo dels of H 11 
region s using the photoionization code Cloudy. iMcGaugbl 
1)19911 ) introduced the y parameter to derive (together with 
the R23 ratio) the oxygen abundance using only the bright 
oxygen lines. The analytical expressi ons for these models 
were given bv lKobulnickv et al.l l| 19991 ). 

It is cl ear from Fig, lc, and as emph asized by 

iDopita et all (|2000l ) and IKewlev fc Dopital (|2002l ). that the 
[O ill] /[On] ratio is probably the best one to use to deter- 
mine the ionization parameter, particularly if the [N n]/[0 11] 
ratio is also measured. To a lesser extent R3 can also be 
used. IKewlev fc Dopital (J2002) provided a procedure to al- 
low a simultaneous solution of both log q and 12 + log(0/H), 
which involves the y and -R23 parameters, although they also 
give a parametrization of the oxygen abundance and log q 
using other parameters such as [Siii]/[Sii] (which is also 
a good estimator of the ionization parameter), N 2O2, N%, 



2004) provided 



S23 and O3./V2. Later. lKobulnickv fc Kewlev 
a parametrization of the IKewlev fc Dopital ( 20021 ) method 
using the y and -R23 paramet ers with a form similar to 
that given by iMcGaughl (1 199 II ) calibration. Following this 
method, log q and the oxygen abundance are computed it- 
eratively using only the bright oxygen lines. 

Pilyugin and his co workers have attempted an em- 
pirical calibration of the R23 and P parameters using the 
bright [Oiii] AA4959,5007 and [On] AA3726,3729 nebular 
lines. To get this calibration, they used a combination of 
photoionization models and empirical data from observed 
H 11 regions where the auroral [O ill] A4363 line was avail- 
able. Indeed, the so-called Pilyugin parameter, P, was in- 



Table 5. Results of the comparison between the oxygen abun- 
dance given by several SEL methods and the oxygen abundance 
assumed by the theoretical models. We indicate the parameters, a 
and 6, of the linear fit xsel = a + ^ x Tnode.l > the correlation coeffi- 
cient, r, of this linear fit, the dispersion of the data, a, and the av- 
erage value of the difference between the abundance given by the 
SEL method and that assumed by the model. In the first row, we 
also include the comparison with the oxygen abundance derived 
following the T c method. The description of the SEL methods is 
the same than that used in Table [4] 



Te 


-0.02 ± 0.14 


1.001 ± 0.017 


0.9959 


0.05 


-0.01 


M91 


-0.20 ± 0.17 


1.005 ± 0.020 


0.9940 


0.07 


-0.15 


KD02 


0.16 ± 0.21 


0.981 ± 0.025 


0.9897 


0.10 


0.00 


KN202 


0.94 ± 0.18 


0.887 ± 0.020 


0.9960 


0.05 


-0.07 


KK04 


0.77 ± 0.10 


0.916 ± 0.012 


0.9975 


0.07 


0.07 


KD 


0.61 ± 0.16 


0.922 ± 0.019 


0.9932 


0.09 


-0.05 


P01 


1.30 ± 0.64 


0.814 ± 0.075 


0.8890 


0.31 


-0.25 


PT05 


1.58 ± 0.51 


0.752 ± 0.061 


0.9113 


0.28 


-0.50 


PVT10 


-1.58 ± 0.24 


1.134 ± 0.029 


0.9901 


0.14 


-0.44 


PM11 


-2.46 ± 0.44 


1.229 ± 0.052 


0.9736 


0.25 


-0.52 


D02 


2.11 ± 0.45 


0.736 ± 0.055 


0.9390 


0.20 


-0.04 


PP04a 


3.43 ± 0.35 


0.575 ± 0.043 


0.9392 


0.25 


-0.04 


PP04b 


3.39 ± 0.43 


0.578 ± 0.052 


0.9141 


0.26 


-0.05 


PP04c 


1.81 ± 1.04 


0.767 ± 0.116 


0.8548 


0.19 


-0.25 


S23 


-1.47 ± 0.66 


1.186 ± 0.081 


0.9481 


0.24 


0.04 


S23023 


-1.76 ± 0.67 


1.208 ± 0.076 


0.9641 


0.17 


0.07 















troduced by IPilvuginl (|2000), after confirming that the 
R23 parameter has a systematic error depending on the 
hardness of the ionizing radiation. Hence, the excitation 
parameter P also considers the ionization degree of the 
Hn region. His first empirical calibrations involv ing both 
the R2 3 and the P paramet er were presented in IPilvuginl 
l|2001ah and IPilvuginl (l2001bl) for high- and low-meta llicity 
Hn regions. Subsequently, IPilvugin fc Thuanl 1)20051 ) and 
IPilvugin. Vilchez fc Thuanl 1 20101 ) tried to improved such 
calibrations including more spectroscopic measurements of 
H 11 regions in spiral and irregular galaxies d erived using the 
T -method. Finally, IPilvugin fc Mattssonl (|201ll ) provided 
another empirical calibration considering only the R3, N2 
and S2 parameters (the so-called NS calibration) , which was 
derived to estimate oxygen abundances in galaxies whose 
spectrum is lacking of the [0 11] AA3726,3729 emission lines. 



4.1 Comparison of the SEL Techniques 

We have applied all of these commonly- used SEL techniques 
to our model spectra to derive oxygen abundances, again 
via a double-blind procedure. The different SEL techniques 
fall into various classes, depending on which combination 
of strong line ratios are used. Table IA3I lists the values 
of all these parameters derived for each model. This table 
also includes the value derived for the q p arameter obtained 
from t he o ptimal calibration provided by IKewlev fc Dopital 
l|2002l ) and iKobulnickv fc Kewlevl (|2004J ). Table H compiles 
the oxygen abundances derived for each model, indicating 
which branch (high or low metallicity) the model spectra 
would fall into when using the .R23 parameter. We considered 
the values of the oxygen abundances derived from the T B 
method to choose the expressions of the lower or the upper 
branch provided by the empirical calibrations. For models 
F, U, W and AF, for which we obtained oxygen abundances 
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Figure 8. Comparison between the total abundanc es used for the models (x-a xis) and those computed using the empirical calibrations 
considering the JV2 and the O3N2 parameters from IDenicolo et al.l 1120021 , and iPettini fc Pagell d2004r ). The continuous line represents 
y = x. Red diamonds and blue triangles represent high- and low-metallicity models, respectively. The vertical dotted line indicates the 
lower (if using the JV2 parameter) or upper (when using the O3N2 parameter) limit of validity of the fits. The dashed line indicates the 
fit to the data within the appropriate metallicity range. 



of 12+log(0/H)~8.4 using the T e method, the results pro- 
vided by both the low- and high- metallicity branches are 
tabulated, and both sets of data are plotted in Figs.[8lto llll 
Note that the calibrations that invoke the JVb, O3N2, S23 or 
■523023 parameters give t he same result for the low and high 
metallicity branches. The lKewlev fc Dopital (|2002l , calibra- 
tion using the N2O2 parameter is only valid for objects in 
the high-metallicity branch. 

To quantify the goodness of these SEL techniques, we 
have performed a linear fit between the oxygen abundances 
derived from these empirical methods and the oxygen abun- 
dances assumed by the models. In all cases, for models with 
12+log(0/H)=8.39 we have assumed the average value be- 
tween the high and low metallicity branches. We note that 
not all models have been considered in this fit, but only 
those within the validity range of each SEL calibration. Ta- 
ble [S] compiles the results for these linear fits, including the 
correlation coefficient, the dispersion of the data and the av- 
erage value of the difference between the abundance given by 
the SEL method and that assumed by the models (offset). 
First row of Table [5] lists the results of a linear fit between 
the oxygen abundances derived following the T e method and 
those imposed by the models. This fit does not consider the 
models with I2+log(0/H)=9.39 (see Fig.©. 



4-1.1 N2 and O3N2 methods 

Let us first consider the techniques which rely upon t he N2 
parameter; that of IDenicolo. Terl evich & Tcrlcvich (2002) 
and those drawn from the lPettini fc Pagell (J2004) which in- 
volve a linear interpolation of the N2 parameter, a cubic fit 
to this same parameter, and their fit using the O3N2 param- 
eter. These are shown in Fig. [5] All these calibrations have 
similar scatter ~ 0.20 — 0.25 dex. However, they are well- 
behaved (if rather curved) in their abunda nce sequence. This 
type o f curvature was also observed by iKewlev fc Ellison! 
(2008), if we can interpret the mass sequence of galaxies as 
equi valent to a n abu ndance sequence. 

lYin et all (120071 ) indicated that the N 2 and the 03^2 
indices are only useful for calibra ting metallicities of galax- 
ies wi th 12+log(0/H)<8.5, while IPerez-Montero fc Continil 
(2009) established that empirical calibrations using 
the O3N2 parameter are not valid for objects with 
12+log(O/H)<8.0. Following Fig. [8] we suggest that the 




.*/ 



12 + log (0/H) [ Model ] 

Figure 9. Comparison between the total abundances used for 
the models (x-axis) and those computed using the empirical cal- 
ibrations considering t h e 523 an d the S23O23 parameters from 
IPerez-Montero fc Diaa J2005I V The continuous line represents 
y = x. Red diamonds and blue triangles represent high- and low- 
metallicity models, respectively. The vertical dotted line indicates 
the upper (using the 523 parameter) or lower (using the S23O23 
parameter) limit of validity of the fits. The dashed line indicates 
the fit to the data within the appropriate metallicity range. 



iV 2 method should not be applied for 12 +log(Q/H)> 8.7- 
9.0, w hile the O3N2 calibration provided bv lPettini fc Page] 
l|2004h is only valid for 12+log(0/H)> 8.7. Importantly, 
in all cases the result has a high uncertainty, ~0.25 dex, 
which is a consequence of lacking of any parameter which 
considers the ionization degree of the gas. Furthermore, in 
the case of galaxies showing an overabundance of nitro- 
gen (e.g. Pustilnik et al. 2004; Lopez-Sanchez et al. 2007; 
Lopez-Sanchez & Esteban 2010; Monreal-Ibcro et al. 2010), 
the application of any N2 or O3N2 calibration will provide 
misleading oxygen abundances. 



4-1.2 S23 and S23O23 methods 

Figure [9] compares the oxygen abundances assumed 
by our models with the results provided using the 
IPerez-Montero fc Diaa (|2005j) calibrations, that consider 
the 523 and S23O23 parameters. In the case of using 
the S23 index, the calibration seems to be valid up to 
12+log(O/H)~9.0, having a dispersion of 0.24 dex. How- 
ever, this method cannot be applied for 12+log(O/H)>9.0. 
On the other hand, the S23O23 calibration works quite well 
at high metallicities, 12+log(0/H)>8.5. But it is not valid 
for 12+log(0/H)<8.2 as the result basically does not de- 



16 A.R. Lopez- Sanchez et al. 




9.5 














/+<'* 


~ 9.0 
M 8.5 












'l 


/ . A : 


2 
















£ 8.0 










1 




: 


5 










fc' 






e 7.5 








A 






- 


™ 7.0 






4/ 


,* 


* 




i 


6.5 




V 


* 


A 









7.5 8.0 8.5 9.0 9.5 

12 + log (O/H) [ Mo, 



7.0 7.5 6.0 8.5 9.0 9.5 

12 + log (0/H) [ Mo, 



6.5 7.0 7.5 8.0 6.5 9.0 9.5 

12 + log (O/H) [ Model ] 



Figure 10. Comparison between the abundances used for the models (x-axis) and those computed via the Pil yugin method us ing 
the R23 and P pa r amete r s. From left to right, pan els show the results for the empirical calibrations provided by iPilvuginl J2001al lbr) , 
iPilvugin fe Thu an (2005). [Pilyugin et al.l J2010T) and lPilvugin fc Ma ttsson (2011J). The continuous line represents y = x. Red diamonds 
and blue triangles represent high- and low-metallicity models, respectively. The dashed line indicates the fit to the data. Note that, in 
the left panel, the three data points with SEL abundances very mu ch larger to thos e of the models are in the intermediate-metallicity 
regime and also have a relatively large logq. Hence, following the IPilvuginl (2001a b) method the high-metallicity calibration should 
be considered and averaged to the value derived using the low-metallicity calibration to get a more appropriate result of the oxygen 
abundance. 
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Figure 11. Comparison b etween the total abunda nces used for the models (x - axis) and those computed using the iMcGaughl 1119911 ) 
method (left panel) and the lKewlev fc Do pita (2002) and iKobulnickv fc Kewlevl 1120041 ) techniques, all of them based on photoionization 
models. The continuous line represents y = x. Red diamonds and blue triangles represent high- and low-metallicity models, respectively. 
The dashed line indicates the fit to the data. The far right panel shows the comparison with the oxygen abundance considering the 
N2O2 method for the h igh-metallicity branch, 12+log(0/H)Js8.35, and the R23 and qKD02o f° r the low-metallicity branch, following 
iKewlev fe Dopital (|2002h . 



pend on metallicity, having a constant value which depends 
on q, and a high dispersion, >0.6 dex. 

Much of the scatte r in the 

iDenicolo. Terlevich fc Terlevichl 120021) . [Pettini fc Pagell 
( 2004 ) and iPerez-Montero fc Diaz] i|2005l ) calibrations is 
due to a failure to include the effects of the ionization 
parameter on the line ratios. That is specially important 
in the analysis of 2D spectroscopic data, as the oxygen 
abundance maps derived using these calibrations may 
directly reflect the ionization structu re of the Hn region 
and n ot actual changes in metallicity (|L6pez-Sanchez et al.l 
120111 ). It is clear that the scatter could be reduced to 
~0.1 dex if the ionization parameter is taken into account, 
since the vertical separation of the points on the figures 
is almost entirely the result of the different ionization 
parameters used in the models. However, this would not 
solve the curvature issue. 



4-1.3 Pilyugin method 

We next consid er the Pi l yugin et al. empirical calibrations 
of R™ and P dPilyuginl feoOlallbl; IPilvugin fc Thuanl 120051 ; 
IPilvugin et alj|20irj|piiyugin fc Mattssonll201ll ). As we dis- 
cussed before, P is an empirical parameter that traces 



the ionization degree of the Hn region. The results we 
obtain are shown in Fig. [10] The main characteristic of 
these SEL calibrations is that there is a systematic off- 
set in the sense that the Pilyugin method tend s to un- 
derest imate the abun dance by up to ^0.2 5 dex l|Pilvuginl 
l2001a|]bT). ~0.4 dex (IPilvugin et all 120101) an d ~0.5 dex 



I Pilvugin fc Thuanl [20051 ; IPilvugin fc Mattssonll201ll) It 

jjbj) ar 



is 
also clear that the scatt er for the IPilvugin (2001a by and 
IPilvugin fc Thuanl (J2005J) empirical calibrations is typically 
0.3 dex. In both calibrations, the scatter becomes consider- 
ably larger in the 12+log(O/H)~8.0-8.4 range, as a result 
of the ambiguity in determining the appropriate branch on 
which the -R23 is located. We note that for this intermediate- 
metallicity regime it is common to give an average value 
between the oxygen abundances derive d for the high a nd 
low metallicity branches if using the IPilvuginl (l2001allbl) 
and Pilyugin fc Thuanl (120051) method (e.g. lMoustakas et al.l 
|2010| ; iLopez-Sanchez fc Estebanll20ich . 

In the case of the IPilvugin et al.l (|2010l ) and 
IPilvugin fc Mattssonl (|201ll ) methods (that also considers 
the N2 and S2 indices), three calibrations are provided for 
low, intermediate and high metallicity. In these cases, the 
scatter is slightly lower, ~0.14 dex and ~0.25 dex, respec- 
tively, than the scatter found using the previous calibrations. 
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The offset is also large, 0.4-0.5 dex, when comparing with 
that provided by the models. 

The very interesting point here is that Pilyugin method 
is giving the best results when it is applied for real ob- 
jects for which the oxygen abund ance has been com- 
puted using the T e me t hod (e .g. iBresolin et ahl 2009; 
lLopez-Sanchez fe Estebanl |2010| ; iMoustakas et all |2010| ). 
However, this method is clearly failing when predicting the 
metallicities of the theoretical models using the bright emis- 
sion lines. We will further discuss this issue in Sect. 14.21 



4-1-4 Kewley and Dopita method 

Although Pilyugin method considers a parameter which is 
related to ionization degree of the gas, the scatter pro- 
vided by this calibration is still relatively high. However, 
the scatter is much reduced when using a SEL method 
which is ba sed on photoion i zation models, as those pro- 
vided by thelMcGaughl (1 1 99 ll ). IKewlev fc Dopital (|2002l ) and 
iKobulnickv fc Kewlevl (J2004J) methods. All these models ex- 
plicitly consider an ionization parameter, y. Fig. [11] shows 
their results. As we see, in all these cases the agreement is 
excellent. We remind the reader that we are plotting both 
the high and low metallicity results for the models with 
12+log(0/H)=8.39. 

Left panel of Fig. [11] shows the case of the McG augh 

(|l99ll ) calibration using the IKobulnickv et all (|l999T ) 
parametrization. This calibration systematically underesti- 
mates the abundances given by the models by ~0.15 dex, but 
their dispersion is only 0.07 dex. In the case of us i ng the 
iKewlev fc Dopital (|2002l ) or IKobulnickv fc Kewlevl (|2004l ) 
methods considering the R23 and y parameter, the results 
are much better. Their method allows to iteratively deter- 
mine the ionization parameter and the oxygen abundance. 
Both, the techniques by Kewley & Dopita and by Kobul- 
nicky & Kewley, eliminate th e scatter in the SEL tech nique 
to a very large extent. For the lKewlev fc Dopital (120021^ tech- 
nique it is ~0.1 dex, while for the IKobulnickv fc Kewlevl 
(2004) technique it is ~0.07 dex excepting the ambiguous 
region around 12+log(0/H)=8.4. 

Far right panel of F ig. 1111 shows a fit u sing t he N2O2 
calibration provided by IKewlev fc Dopital ([2002]) for the 
high- metallicity branch, 12+log(0/H)^8.35, and relied on 
the R23 and the technique they use for obtaining the ioniza- 
tion parameter for the lower metallicities. This assumption 
provides very good results, as this method has a very low 
offset, 0.05 dex, and dispersion, 0.09 dex. 

The slight offset in t he ab undances derived follow i ng the 
IKewlev fc Dopital (|2002l ) and IKobulnickv fc Kewlevl (|2004l ) 
methods is probably due to the fact that they used a differ- 
ent grid of models. The mode ls we use in this paper are more 
closely related to the earlier IKewlev fc Dopital ([2002]) grid, 
except in respect of the much more sophisticated treatment 
of dust physics employed here. 

However, the recent analysis of the ionized gas within 
a sample of ^40 strong sta r -form ing galaxies performed 
by lLopez-Sanchez fc Estebanl (J2010f) shows that t he oxygen 
abundances provided by phot o ionization models ( McGaugh 
Il99ll ; IKewlev fc Dopital |2002| ; IKobulnickv fc Kewlevl I2OO4T ) 
are, systematically, 0.2-0.3 dex higher than the oxygen 
abundances provided by the T e method. Th i s observational 
result has been also found bv lYin et all (J2007T ). IBresolin et al.l 
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Figure 12. Comparison between the oxygen abundance derived 
following the T c , RL, P and KD methods for deep, high-resolution 
observations of Galactic (open symbols) and extragalactic (filled 
symbols) Hn regions for which oxygen recombination lines have 
been detected. The ai-axis plots the T e abundance, where the y- 
axis represents the difference between RL (black squares), P (blue 
diamonds) or KD (red triangles) abundances and the T e abun- 
dance. The black dotted line plots the average difference between 
the T and the RL methods. The blue dashed line represents the 
average difference between the T e and the P methods. The red 
dotted-dashed line represents the average difference between the 
T c and the KD methods. 



(2009), IMoustakas et alJ <|2010h and lLopez-Sanchez et all 
(|2011| l : in some objects the differences reach up to 0.6 dex. 
Interestin g, however , the p hotoionization models devel- 
oped by iDors et al.l l|201lT ) provided O/H values close 
to those obtained using the T c method. Detailed tailor- 
made photoionization mo dels of individual galaxies (e.g. 
iPerez-Montero et alj|20ich also provided results similar to 
those derived from the observed electron temperatures of 
the ionized gas. 



4.2 Comparison of the SEL techniques 
with real data 

Hence, we have further investigated the validity of the 
SEL techniques using real data. For this, we aim to use 
the deepest, highest quality observational data of Galac- 
tic and extragalactic Hn regions nowadays available. Par- 
ticularly, we also want to compare the results provided 
by the T e and the most common SEL techniques (the 
Pilyugin and the Kewley & Dopita methods) with the 
oxygen abundances derived using oxyge n recombination 
lines (RL). As it is we l l known (e.g. | Peimbertl 12003 



Esteban et ail |2002|. liool. 120091; IGarcia-Roias et al 



2004, 



20051 ; IGarcia-Roias fc Estebarj|2007l ; iPeimbert et alj|2007t ). 

in all cases the oxygen abundances determined using RL are 
0.2-0.3 dex higher than the oxygen abundances derived us- 
ing the T e method, which is based on collisionally excited 
lines (CEL). The effect does not depend on metallicity, elec- 
tron density, electron and Balmer temperatures or ionization 
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Table 6. Results of the oxygen abundances given by the Te, RL, P and KD methods for Galactic (G) and extragalactic (E) Hn regions 
for which oxygen recombination lines have been detected. 



Object 



Rs 



logiV 2 



T e 



Typo 



M42 

NGC 3576 
S311 
M16 
M20 

NGC 3603 
M8 
M17 

30 Dor 

NUB 

N66 

NGC 6822V 

NGC5253A 

NGC5253B 

NGC5253C 

NGC 595 

NGC 604 

H1013 

NGC 5461 

VS44 

NGC 2363 

K932 



2.981 
4.169 
1.318 
0.668 
0.350 
0.828 
3.038 
3.038 

1.823 
1.925 
2.052 
3.000 
2.429 
1.538 
0.554 
0.596 
0.255 
1.038 
0.121 
1.603 
4.134 
1.085 



0.375 
2.447 
5.591 
3.547 
0.009 
0.389 
5.002 
5.002 

1.109 
5.728 
6.349 
3.764 
3.860 
0.091 
1.229 
0.287 
0.032 
1.178 
0.030 
0.690 
1.111 
2.203 



-1.247 
-0.767 
-1.720 
-2.055 
-0.549 
-1.932 
-1.184 
-1.184 

-0.692 
-1.502 
-1.360 
-1.296 
-1.020 
-0.430 
-2.160 
-0.792 
-0.627 
-0.938 
-0.865 
-1.618 
-1.096 
-1.842 



8.51±0.03 
8.56±0.03 
8.39±0.05 
8.56±0.07 
8.53±0.06 
8.46±0.05 
8.51±0.05 
8.52±0.04 

8.33±0.02 

8.41±0.00 
8.11±0.00 
8.08±0.03 
8.18±0.04 
8.19±0.04 
8.28±0.04 
8.45±0.03 
8.38±0.02 
8.45±0.06 
8.41±0.03 
8.36±0.02 
7.76±0.02 
8.41±0.02 



!.71±0.03 

!.74±0.06 
!.54±0.10 
!.81±0.07 
!.71±0.07 
!.72±0.05 
!.71±0.04 
!.76±0.04 

!.54±0.05 
!.74±0.00 
!.47±0.00 
!.37±0.09 
!.42±0.13 
!.37±0.10 
!.53±0.09 
!.69±0.05 
!.60±0.03 
!.73±0.09 
!.49±0.06 
!.61±0.04 
!.04±0.05 
!.62±0.03 



8.57 
8.51 
8.28 
8.45 
8.45 
8.45 
8.50 
8.58 

8.39 
8.43 
8.10 
8.09 
8.12 
8.12 
8.30 
8.51 
8.61 
8.62 
8.56 
8.59 
7.71 
8.52 



8.77 
8.80 
8.65 
8.79 
8.93 
8.68 
8.95 
8.78 

8.63 
8.80 

8.40 
8.39 
8.48 
8.47 
8.59 
8.72 
8.73 
8.90 
8.73 
8.72 
8.10 
8.77 



G 
G 
G 
G 
G 
G 

c; 

G 

E 
E 
E 
E 
E 
E 
E 
E 
E 
E 
E 
E 
E 
E 



E04 
GR04 
GR05 
GR06 
GR06 
GR06 
GR07 
GR07 

Pe03 
T03 
T03 
PPR05 
LS07 
LS07 
LS07 
E09 
E09 
E09 
E09 
E09 
E09 
E09 



References: E04: | Esteban et all (20041), E09:lEsteban et al.l i 20091 ). GR04: Garcia-Rojas et al. (2004(1 GR05: jGarcia-Rojas et al. 12005), GR06: 
Gareia-Roias et al .1 (20060, g"r07 : Garcia-Rojas & Esteban 12003), LS07: iLdpez-Sanchez et al.l 720071 1. Pe03: jPcimbert 12003), PPR05: 



Peimbert, Peimbert fc RuiJ (20051) , T03: iTsamis et al.U 2 



degree l|Garcia-Roias fc Estebanl 120071 ; iMesa-Delgado et alj 

2009a) and it may be attributed to the presence of temper- 
ature fluctuations in Hn regions. 

Moreover, r ecent analyses (Peimb ert. Peimbert fc Ruia 



2005: 



Przvbilla. Nieva fc Butlerl 



200c 



Simon-Diaz fc Stasiriskal 1201 ll ) have found that the 



oxygen nebular abundance based on RL agrees better with 
the abundances of the stars associated with the nebulae 
than the one derived from the T c method, suggesting that 
RL abundances are more reliable than the abundances 
derived from CEL. 

Table [6] lists all the 22 Galactic and extragalactic H II 
regions found in the literature for which oxygen abundances 
are available using both the RL and the T c methods. In- 
deed, the average offset between T c and RL abundances 
is 0.22 dex. Table [S] also compiles the derived -R23, P , y 
and JV2 indexes derived from the reddening-corrected data. 
We must emphasize that, in almost all the cases, the er- 
ror in the [On] AA3726,3729 lines is less than 5%. Using 
these observational parameters, we then applied the Pilyu- 
gin et al. method and the Kewley & Dopita method to de- 
termine the oxygen abundances using these SEL te chniques- 
Table [6 ] co mpiles the average value d erived for the IPilvuginl 
|2001al lbh: IPilvugin fc Thuanl |2005h ; iPilvugin et all (|201ff ) 
techniques (P) and th e average value obtained for the 
iKewlev fc Dopital j|2002h (using the .R 2 3 , y a nd N2O2 param- 
eters) and the lKobulnickv fc Kewlevl f200_4) method (KD). 



As we see in Table [6l we find again that the results 
obtained using the P technique agree quite well with the 
oxygen abundance determined using the T e method (offset 
of 0.03 dex), but the results provided by the KD method 
are systematically ~0.30 dex larger. That is exactly the op- 
posite behavior we discussed in the previous section. This 
result is graphically shown in Fig. 1121 which compares the 
Te abundance (x-axis) with the oxygen abundances derived 



from RL (black squares), the P method (blue diamonds) and 
the KD technique (red triangles). 

Interestingly, the results provided using the KD tech- 
nique have a better agreement with the values determined 
from RL. These are systematically ~0.08 dex lower than the 
oxygen abundances given by the KD method, but both agree 
within the errors (the dispersion of the RL data is ~0.06 dex, 
while the dispersion of the KD data is ~0.07 dex). 
Thi s observational result , whi ch was previously noted 
by Peimbert et a.1.1 (2007) and ILopez-Sanchez fc Estebanl 
(2010), is actually suggesting that the real metallicities of 
the H 11 regions may be not those given by the standard 
T e method based on line ratios of collisionally excited lines 
(CELs), but the values derived using recombination lines 
(RL). If this is correct, then the KD method is the only 
SEL technique that is reproducing those metallicities. 

As we have seen along this paper, the oxygen abun- 
dances of our model spectra are well reproduced using 
the T e method and SEL t echni q ues based on photoion- 
ization models jMcGaughl Il99ll ; IKewlev fc Dopital |2002| ; 
iKobulnickv fc Kewlevl 120041 ) but 0.3-0.5 dex higher than 
those estimated using the Pilyugin method. However, the 
situation is inverted in other works based on the compari- 
son of observations w it h the results of t he SEL methods (e.g . 
Peimbert et al.1 120071; lYin et all I2OO7I; iBresolin et al.l 1 2(305 : 



Lopez- Sanchez fc Estebanl 201CI; i Moustak as et al 



2010; 



Lopez-Sanchez et al.ll201ll ; lRosales-Ortega et al.ll201lT ). i.e 
the oxygen abundance in real nebulae is better reproduced 
by those SEL nrethods which rely in an empirical calibration 
derived using H 11 regions where a direct estimation of the 
electron temperature exists. 

A possible explanation of this puzzling situation is that 
current photoionization models are not properly reproduc- 
ing the behavior of collisionally excited lines (CELs) in 
real nebulae. As we said before, several authors have in- 
voked the presence of temperature fluctuations in the neb- 
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ulae (|Peimbertlll967h in order to address the differences be- 
tween the abundances determined from CELs and RLs. Due 
to the strong dependence on T c of the intensity of CELs, 
such temperature fluctuations cause the metallicities de- 
rived from CELs to be underestimated. Although our mod- 
els consider global temperature gradients due to the hard- 
ening of the ionizing radiation field along the nebula, they 
did not reproduce the values of the temperature fluctua- 
tions parameter (£2) that are necessary to make the O/H 
ratios determined from RLs and CELs to converge (e.g. 
iGarcfa-Roias fc Estebanl 120071 ) . These hypothetical fluctu- 
ations must be of small spatial scale and may be produced 
by, for example, density inhomogeneities or localized depo- 
sition of mechanical energy that heat the gas. 

Summarizing, assuming the validity of the scenario out- 
lined above, the oxygen abundances determined from RLs 
and using SEL techniques based on photoionization mod- 
els (KD method) would be closer to the real ones, while 
values determined from CELs and empirical calibrations (P 
method) would provide abundances about 0.2-0.3 dex lower. 
This should have a profound impact in current galaxy metal- 
licity estimations that have been obtained using CELs or 
any empiric al calibration-such as those provided by Pilyu- 
gin et al. or lPettini fc Page] (|2004T ) - that should be revised 
upwards. 

Finally, we want to note that the small offset of 
~0.08 dex between the RL and the KD may be a con- 
sequence of the depletion of the oxygen atoms in dust 
grains. Our models are taking into account this effect, 
but observations suggest that oxygen depletion in dust 
grains may be decreasing the actual gas-phase metallic- 
ity (i.e., that derived from RL) of Hn regions by 0.08- 
0.12 dex llMesa-Delgado et ai1|2009bl ; IPeimbert fc Peimbertl 
l2010l ; ISim6n-Di'az fc Stas iiiska 20 13). 



5 CONCLUSIONS 

In this paper, we have subjected the commonly-used 
techniques for deriving chemical abundances to a critical 
"double-blind" analysis. We find that the originally-used T 
method is capable of delivering reliable abundances for O, 
N, S, and CI. However, the abundances inferred for the no- 
ble gases Ne and Ar are considerably less secure, although 
probably for different physical reasons. 

The electron temperature itself is a reasonably good in- 
dicator of the oxygen abundance, as the one is the dominant 
controller of the other. Error arises because there is also a 
weaker dependence of electron temperature on the ionization 
parameter. For the [O ill] lines, this renders an abundance 
obtained from the electron temperature alone unreliable be- 
low 12+log(0/H) < 8.0. 

The temperature derived from the temperature- 
sensitive collisionally excited line ratios is always somewhat 
different from emission-weighted temperature defined for the 
principal emission lines -A4959+A5007 in the case of [O ill]-. 
This is due to the variability of electron temperature through 
the zone containing the ion of interest due to hardening of 
the radiation field and spherical divergence of the radia- 
tion field (which changes the photoelectric heating rate), 
and change in ionization state (which influences the cooling 
rate of the plasma). In this sense, the models naturally pro- 



duce electron temper ature gradients in s ide H 11 regions, such 
as those discussed in Istasiriskal (|1978T ) , ICarnettJ (|1992T ) and 



IStashiskal (120051 ). 

Although our models can reproduce these "global" 
temperature fluctuations, they cannot simulate the small- 
sca le temperature fl uctua tions, such as those f ir st dis cussed 
by IPeimbert! (|l967n and IPeimbert fc Costerol (|l969h . pro- 
duced by turbulent or photo-evaporative flows, shocks, or 
radiation shielding by dense inclusions or globules of un- 
ionized material. Thus we can be reasonably sure that the 
true temperature fluctuations in real H 11 regions should be 
greater than the models predict. The effect of this would 
be to raise the temperature as measured by the [O ill] 
(A4959+A5007)/A4363, which would then feed into a sys- 
tematic underestimate in the oxygen abundance as deliv- 
ered by the T c method. Hence, accepting the existence of 
small-scale temperature fluctuations within the ionized gas 
and assuming that the RLs give the true metallicities, those 
oxygen abundances determined from CELs and empirical 
calibrations (P method) would provide abundances which 
are 0.2-0.3 dex lower than the real ones. That is because of 
the strong dependence of the intensity of CELs on T e : the 
existence of such small-scale temperature fluctuations will 
cause an underestimation of the oxygen abundances using 
CELs lines or strong-line methods which are based on the 
bright nebular lines. However, oxygen abundances derived 
using SEL techniques based on photoionization models (such 
as the KD method, for which the T c is fixed a priori) would 
be closer to the real values. 

We have demonstrated that those SEL techniques 
which do not explicitly solve for the ionization param- 
eter are quite unreliable, resulting in large scatter, sys- 
temic error, and abundance-dependent errors. To obtain 
the ionization parameter, we need a measurement of 
the [Oni] A5007/[On] AA3726,3729 ratio. For abundances 
12+log(0/H)< 8.0 the oxygen abundance could then be es- 
timated from the [O ill] A5007/H/3 line ratio. This is verified 
in Fig 1 131 which provides a new diagnostic plot for this low 
abundance regime. 

At high abundance, the most sensitive diagnos- 
tic plot is that of [Oni] A5007/[On] A3726, 3729 vs. 
[Nil] A6584/[Oll] AA3726.3729, shown in panel c of Fig. Q] 
The main error here is the degree to which the N/O abun- 
dance has been accurately calibrated against the O abun- 
dance. We believe that the relationship given in Eq. [T] 
is probably good, sin c e the model grids reproduce the 
IVeilleux fc Osterbrockl (119871 ) diagnostic plots shown in 
Fig. [1] panels (a) and (b). 

Figures [TJ panel (c) and [T3] thus provide a SEL tech- 
nique that is applicable to the determination of both metal- 
licity and ionization parameter over the full abundance 
range, with the caveat that the sensitivity of this technique 
is rather poor in the range 8.0 < 12 + log(0/H) < 8.4. In 
this range, common to all SEL techniques, rather accurately- 
measured line ratios are required. 
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APPENDIX A: ELECTRON TEMPERATURES, 
ELECTRON DENSITIES AND IONIC 
ABUNDANCES FOR THE MODEL H n 
REGIONS. 

Here we compile the derived parameters for the model H II 
regions. The first table, Table |A"T1 gives the derived densities 
and temperatures in the various zones of the H II region. Ta- 
ble IA2I gives the derived ionic abundances. Finally, Tabic 
IA3I gives the line ratio values used to compute the oxy- 
gen abundance using the empirical calibrations as well as 
the ionization paramete rs der ived for the models us i ng the 
iKewlev fc Dopital (J2002T I and iKobulnickv fc Kewlevl (J2004T I 
techniques. 
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Table Al. Electron density and temperature obtained from the analysis of emission line ratios of model spectra using IRAF task ternden. 









n e [ 


cm 2 ] 






T 


High [ K ] 






T c 


Low [ K 




Model 


[Sn] 


[On] 


[Ni] 


[CI III] 


Adopted 


[Oin] 


[Sin] 


[Ar in] 


Adopted 


[On] 


[Sn] 


[Nil] 


Adopted 


A 


6 


21.3 


12 


113 


86 


8465 


8037 


8389 


8297 


8419 


6419 


9274 


8037 


B 





213 


7 


136 


96 


14947 


12222 


13307 


13492 


12141 


7621 


13658 


11140 


C 


.30 








68 


49 





4425 


4519 


4472 





3843 


4680 


4261 a 


D 





236 


7 


135 


126 


12663 


11009 


11706 


11792 


11122 


7337 


12820 


10426 


E 





228 


8 


130 


122 


11459 


10134 


10732 


10775 


10311 


7086 


11802 


9733 


F 





226 


9 


130 


121 


10498 


9625 


10159 


10094 


9834 


7023 


11254 


9370 


G 


25 


182 


2.3 


76 


76 


4855 


4993 


4908 


4918 


5423 


4535 


5660 


5206 


fi 


.33 








68 


50 


4285 


4210 


4236 


4243 





3737 


4474 


4105 a 


1 


19 


196 


21 


87 


80 


5830 


5873 


5973 


5892 


6424 


5313 


6986 


6241 


.1 


35 








71 


53 





3910 


3768 


3839 





3633 


4269 


3951 a 


K 





231 


7 


1.35 


125 


12413 


10843 


11519 


11591 


10968 


7304 


12633 


10301 


L 





229 


8 


130 


122 


11394 


10139 


10737 


10756 


10297 


7122 


11796 


9738 


M 





210 


7 


138 


128 


13637 


11749 


12613 


12666 


11614 


7514 


13408 


10845 


N 


6 


211 


12 


113 


86 


8129 


7811 


8144 


8028 


8409 


6435 


9280 


8041 





21 


182 


21 


77 


75 


5769 


5542 


5665 


5658 


5681 


4660 


5996 


5445 


P 





211 


7 


138 


128 


14096 


12034 


12953 


13027 


11833 


7556 


13658 


11015 


Q 





239 


7 


135 


127 


13718 


11640 


12568 


12642 


11462 


7454 


13202 


10706 


R 


16 


191 


19 


86 


78 


6225 


6158 


6283 


6222 


6549 


5322 


6996 


6289 


S 





232 


8 


135 


125 


11890 


10586 


11248 


11241 


10661 


7284 


12268 


10071 


T 





231 


8 


136 


126 


12181 


10781 


11506 


11489 


10804 


7314 


12451 


10189 


U 





228 


9 


131 


122 


10525 


9641 


10182 


10116 


9913 


7044 


11364 


9440 


V 





211 


7 


139 


129 


14562 


12172 


13223 


13319 


11840 


7525 


13640 


11001 


w 





225 


9 


129 


121 


10537 


9589 


10116 


10080 


9748 


6962 


11119 


9276 


X 


10 








78 


59 





3930 





3930 





3695 


4338 


4016 a 


Y 





213 


7 


136 


128 


15740 


12547 


13728 


14005 


12101 


7554 


13892 


11182 


Z 


13 


191 


17 


90 


78 


6941 


6580 


6790 


6770 


6768 


5356 


7252 


6458 


AA 


22 


181 


2.3 


75 


75 


5364 


5295 


5338 


5332 


5498 


4588 


5788 


5291 


AB 





231 


8 


1.35 


125 


12452 


10927 


11628 


11669 


11001 


7359 


12678 


10346 


AC 


15 


191 


17 


89 


78 


6741 


6483 


6654 


6626 


6712 


5378 


7189 


6426 


AD 





243 


6 


1.38 


129 


15240 


12497 


13644 


13793 


12078 


7563 


13888 


11176 


AE 





238 


7 


1.39 


128 


13248 


11484 


12330 


12354 


11363 


7459 


13115 


10645 


AF 





222 


10 


125 


119 


10289 


9308 


9814 


9803 


9470 


6799 


10740 


9003 


AG 


5 


212 


12 


113 


85 


8765 


8182 


8545 


8497 


8477 


6389 


9318 


8061 



a ionic task does not allow electron temperatures lower than 5,000 K for the low ionization ions, so we adopted T c (low) = 5,000 K for 
these cases. 
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Table A2. Ionic abundances derived for the results of the photoionization models using IRAF task ionic considering the electron 
temperature derived from all the high-ionization species, T c (high), and the low ionization species, T c (low), separately. All numbers are 
expressed in the form of 12 + log (X/H). The label disp. indicates the dispersion of the data when considering the values given by 
different emission lines. 



Model 




O + /H+ 


O++/H+ 


N+/H+ 


S+/H+ 


S++/H+ 


Nc+ + /H+ 


Ar+/H+ 


Ar+ + /H + 


C1++/H+ 


Cl 3 +/H+ 


A 




8.717 


8.126 


7.333 


6.327 


7.074 


7.683 


6.458 


4.364 


4.272 


3.578 




disp. 


0.0S1 


0.016 


0.085 


0.095 


0.0JS 


0.07J 


0.007 


0.001 


0.001 


0.0.90 


B 




7.283 


7.609 


5.635 


4.928 


5.873 


6.941 


5.326 


4.915 


2.810 


3.251 




disp. 


0.036 


0.049 


0.084 


0.122 


0.036 


0.071 


0.006 


0.000 


0.001 


0.090 


C 




8.873 


7.101 


8.257 


7.224 


7.365 


7.347 


6.586 


0.000 


4.575 


2.487 




dis-p. 


0.115 


4.100 


0.070 


0.189 


0.011 


0.071 


0.013 


0.000 


0.001 


0.090 


D 




7.489 


6.809 


5.840 


5.208 


5.727 


6.468 


5.132 


2.969 


2.882 


2.502 




disp. 


0.029 


0.039 


0.091 


0.119 


0.028 


0.071 


0.004 


0.001 


0.000 


0.0.90 


E 




8.173 


6.915 


6.606 


6.003 


6.282 


6.884 


5.703 


2.871 


3.517 


2.695 




disp. 


0.029 


0.0,97 


0.092 


0.113 


0.02S 


0.071 


0.00,9 


0.001 


0.000 


0.090 


F 




8.252 


8.127 


6.717 


5.796 


6.752 


7.544 


6.157 


4.787 


3.868 


3.666 




disp. 


0.026 


0.026 


0.091 


0.106 


0.625 


0.071 


0.00,9 


0.001 


0.000 


0.090 


G 




9.089 


8.315 


8.038 


6.764 


7.525 


7.803 


6.747 


4.544 


4.725 


3.535 




disp. 


0.0S8 


0.018 


0.079 


0.088 


0.0J5 


0.071 


0.003 


0.000 


0.000 


0.090 


II 




8.710 


7.821 


8.152 


7.003 


7.564 


7.658 


6.802 


3.624 


4.771 


3.019 




disp. 


0.144 


0.017 


0.122 


0.219 


0.008 


0.07J 


0.002 


0.000 


0.00 J 


0.0.90 


I 




8.963 


8.309 


7.709 


6.447 


7.240 


7.754 


6.559 


4.909 


4.416 


3.763 




disp. 


0.061 


0.076 


0.137 


0.058 


0.0-jS 


0.071 


0.041 


0.000 


0.000 


0.090 


J 




8.372 


8.105 


7.892 


6.667 


7.676 


7.710 


6.781 


4.095 


4.855 


3.231 




di,sp. 


0.170 


4. 680 


0.178 


0.238 


0.022 


0.071 


0.028 


0.002 


0.001 


0.0.90 


K 




7.778 


7.081 


6.155 


5.502 


6.028 


6.754 


5.433 


3.223 


3.194 


2.765 




disp. 


0.026 


0.039 


0.087 


0.120 


0.028 


0.071 


0.004 


0.001 


0.000 


0.0.90 


L 




8.313 


7.388 


6.758 


6.015 


6.505 


7.129 


5.905 


3.456 


3.711 


3.024 




disp. 


0.028 


0.035 


0.091 


0.112 


0.027 


0.071 


0.002 


0.001 


0.000 


0.090 


M 




7.536 


7.470 


5.884 


5.167 


6.022 


6.876 


5.431 


4.191 


3.036 


3.096 




disp. 


0.029 


0.038 


0.089 


0.120 


0.029 


0.07J 


0.003 


0.001 


0.000 


0.090 


N 




8.481 


8.378 


7.079 


6.029 


7.065 


7.804 


6.458 


5.041 


4.216 


3.858 




disp. 


0.031 


0.010 


0.085 


0.094 


0.0J6 


0.071 


0.010 


0.001 


0.001 


0.090 


C) 




9.185 


7.426 


8.201 


7.180 


7.307 


7.405 


6.597 


3.204 


4.517 


2.846 




disp. 


0.025 


0.000 


0.047 


0.119 


0.0,9^ 


0.071 


0.018 


0.000 


0.000 


0.0.90 


P 




7.218 


7.183 


5.548 


4.862 


5.712 


6.581 


5.124 


3.946 


2.705 


2.818 




disp. 


0.030 


0.039 


0.089 


0.122 


0.030 


0.071 


0.00,9 


0.001 


0.000 


0.090 


Q 




7.644 


7.880 


6.022 


5.232 


6.224 


7.223 


5.667 


5.094 


3.210 


3.505 




disp. 


0.029 


0.042 


0.090 


0.120 


0.032 


0.071 


0.003 


0.001 


0.000 


0.090 


R 




8.984 


8.338 


7.788 


6.615 


7.368 


7.862 


6.711 


4.579 


4.579 


3.652 




disp. 


0.034 


0.001 


0.082 


0.089 


0.00s 


0.07J 


0.009 


0.000 


0.000 


0.090 


S 




8.110 


7.855 


6.516 


5.674 


6.535 


7.300 


5.927 


4.379 


3.643 


3.414 




disp. 


0.027 


0.0,9,9 


0.090 


0.112 


0.026 


0.071 


0.001 


0.001 


0.000 


0.0.90 


T 




7.896 


8.008 


6.297 


5.444 


6.472 


7.373 


5.896 


5.017 


3.510 


3.603 




disp. 


0.027 


0.033 


0.090 


0.114 


0.027 


0.071 


0.001 


0.00 J 


0.000 


0.0.90 


U 




8.067 


8.227 


6.531 


5.607 


6.692 


7.596 


6.117 


5.290 


3.755 


3.812 




disp. 


0.026 


0.026 


0.091 


0.107 


0.02S 


0.071 


0.00,9 


0.002 


0.000 


0.0.90 


V 




7.349 


7.579 


5.698 


4.983 


5.921 


6.924 


5.369 


4.758 


2.868 


3.221 




disp. 


0.030 


0.043 


0.089 


0.123 


0.0,9,9 


0.071 


0.004 


0.00 J 


0.000 


0.090 


w 




8.481 


7.804 


6.972 


6.106 


6.759 


7.406 


6.150 


3.992 


3.946 


3.351 




disp. 


0.026 


0.029 


0.091 


0.107 


0.02-j 


0.071 


0.002 


0.00 J 


0.000 


0.090 


X 




8.199 


7.632 


7.698 


6.368 


7.457 


7.095 


6.138 


3.674 


4.599 


2.834 




disp. 


0.152 


4-407 


0.159 


0.227 


0.000 


0.071 


0.607 


0.003 


0.000 


0.090 


Y 




7.004 


7.312 


5.337 


4.657 


5.561 


6.641 


5.015 


4.628 


2.484 


2.956 




disp. 


0.031 


0.05,9 


0.088 


0.126 


0.039 


0.071 


0.009 


0.00 J 


0.000 


0.090 


z 




9.179 


7.161 


8.024 


7.151 


7.158 


7.350 


6.502 


2.898 


4.387 


2.758 




disp. 


0.036 


0.025 


0.087 


0.098 


0.020 


0.071 


0.002 


0.00 J 


0.00 J 


0.090 


AA 




9.292 


8.103 


8.240 


7.056 


7.484 


7.746 


6.769 


4.105 


4.698 


3.336 




disp. 


0.0,97 


0.008 


0.083 


0.090 


0.006 


0.071 


0.001 


0.000 


0.000 


0.090 


AB 




8.016 


7.570 


6.405 


5.639 


6.349 


7.093 


5.742 


3.905 


3.478 


3.168 




disp. 


0.028 


0.0,96 


0.090 


0.116 


0.027 


0.071 


0.002 


0.001 


0.000 


0.090 


AC- 




9.138 


7.906 


7.964 


6.903 


7.303 


7.666 


6.654 


3.869 


4.529 


3.294 




disp. 


0.035 


0.017 


0.084 


0.093 


0.0J6 


0.071 


0.00,9 


0.00 J 


0.001 


0.0.90 


AD 




7.042 


7.285 


5.374 


4.690 


5.604 


6.625 


5.055 


4.492 


2.535 


2.931 




disp. 


0.031 


0.047 


0.089 


0.125 


0.035 


0.071 


0.005 


0.00 J 


0.000 


0.090 


AE 




7.766 


7.814 


6.141 


5.348 


6.304 


7.189 


5.722 


4.717 


3.322 


3.422 




disp. 


0.028 


0.0,97 


0.090 


0.118 


0.029 


0.071 


0.002 


0.00 J 


0.000 


0.0.90 


AF 




8.585 


7.112 


7.118 


6.439 


6.669 


7.157 


6.084 


3.024 


3.914 


2.857 




disp. 


0.027 


0.033 


0.091 


0.107 


0.026 


0.071 


0.001 


0.00 J 


0.000 


0.090 


AG 




8.851 


7.570 


7.503 


6.639 


7.005 


7.456 


6.400 


3.515 


4.239 


3.153 




disp. 


0.033 


0.024 


0.085 


0.098 


0.021 


0.071 


0.004 


0.001 


0.001 


0.090 
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Table A3. Parameters used to compute the oxygen abundance using empirical calibrations. 



R-2 



R 3 



f-2:i 



logiV 2 lo^OsN 2 log02iV 2 



l°g9J<rD02° 



1OS1KK04 



A 


4.169 


2.447 


6.616 


0.370 


-0.231 


-0.640 


1.026 


1.387 


0.623 


2.445 


7.72 


7.32 


B 


0.668 


3.547 


4.215 


0.842 


0.725 


-1.928 


2.481 


1.880 


0.060 


0.367 


8.40 


7.87 


C 


0.350 


0.009 


0.359 


0.025 


-1.583 


-0.422 


-1.619 


0.093 


1.001 


1.560 


6.69 


6.47 


D 


0.828 


0.389 


1.216 


0.319 


-0.329 


-1.805 


1.397 


1.850 


0.098 


0.269 


7.16 


7.13 


E 


2.981 


0.375 


3.356 


0.112 


-0.901 


-1.120 


0.691 


1.721 


0.510 


1.028 


7.30 


6.87 


F'' 


3.038 


5.002 


8.040 


0.622 


0.217 


-1.057 


1.754 


1.667 


0.284 


1.960 


7.81 


7.65 , 7.74 


G 


0.596 


0.287 


0.883 


0.325 


-0.318 


-0.665 


0.120 


0.567 


0.352 


1.452 


7.47 


7.83 


II 


0.255 


0.032 


0.288 


0.112 


-0.898 


-0.500 


-0.993 


0.034 


0.620 


1.337 


6.78 


7.52 


I 


1.038 


1.178 


2.215 


0.532 


0.055 


-0.811 


0.879 


0.954 


0.266 


1.604 


7.78 


8.09 


J 


0.121 


0.030 


0.152 


0.200 


-0.603 


-0.738 


-0.783 


-0.051 


0.294 


0.875 


6.94 


8.58 


K 


1.603 


0.690 


2.293 


0.301 


-0.366 


-1.491 


1.328 


1.823 


0.191 


0.522 


7.08 


7.14 


L 


4.134 


1.111 


5.245 


0.212 


-0.571 


-0.969 


1.013 


1.712 


0.524 


1.384 


7.61 


7.09 


M 


1.085 


2.203 


3.288 


0.670 


0.308 


-1.715 


2.056 


1.877 


0.098 


0.481 


8.00 


7.56 


N 


2.429 


3.860 


6.289 


0.614 


0.201 


-0.893 


1.477 


1.405 


0.313 


1.959 


7.90 


7.85 


O 


1.538 


0.091 


1.629 


0.056 


-1.227 


-0.303 


-0.739 


0.617 


1.397 


2.610 


6.87 


6.99 


P 


0.554 


1.229 


1.782 


0.689 


0.346 


-2.033 


2.120 


1.903 


0.050 


0.247 


7.54 


7.54 


Q 


1.318 


5.591 


6.909 


0.809 


0.627 


-1.593 


2.338 


1.840 


0.110 


0.721 


8.31 


7.90 


R 


1.823 


1.109 


2.932 


0.378 


-0.216 


-0.565 


0.608 


0.953 


0.532 


2.210 


7.55 


7.73 


S 


3.000 


3.764 


6.764 


0.556 


0.099 


-1.169 


1.742 


1.773 


0.260 


1.262 


7.94 


7.53 


T 


1.925 


5.728 


7.653 


0.748 


0.473 


-1.375 


2.131 


1.787 


0.158 


1.061 


8.34 


7.81 


U c 


2.052 


6.349 


8.401 


0.756 


0.490 


-1.233 


2.034 


1.672 


0.188 


1.356 


8.37 


7.85 , 7.97 


V 


0.744 


3.232 


3.976 


0.813 


0.638 


-1.883 


2.390 


1.882 


0.066 


0.400 


8.32 


7.80 


w c 


4.906 


2.354 


7.260 


0.324 


-0.319 


-0.814 


1.183 


1.632 


0.565 


1.919 


7.62 


7.28 , 7.34 


X 


0.079 


0.012 


0.092 


0.134 


-0.812 


-0.961 


-1.002 


-0.012 


0.145 


0.550 


7.12 




Y 


0.357 


1.968 


2.325 


0.846 


0.741 


-2.242 


2.533 


1.922 


0.033 


0.192 


8.42 


7.82 


z 


3.360 


0.107 


3.467 


0.031 


-1.495 


-0.288 


-0.683 


0.941 


2.049 


3.382 


6.90 


6.76 


AA 


1.085 


0.275 


1.361 


0.202 


-0.595 


-0.434 


-0.128 


0.597 


0.742 


2.104 


7.01 


7.49 


AB 


2.707 


2.180 


4.887 


0.446 


-0.094 


-1.257 


1.593 


1.816 


0.257 


0.959 


7.91 


7.36 


AC 


2.975 


0.545 


3.520 


0.155 


-0.737 


-0.354 


0.088 


0.954 


1.127 


2.874 


7.44 


7.22 


AD 


0.389 


1.796 


2.185 


0.822 


0.664 


-2.206 


2.458 


1.923 


0.035 


0.206 


8.34 


7.76 


AE 


1.710 


4.535 


6.245 


0.726 


0.424 


-1.486 


2.140 


1.846 


0.141 


0.844 


8.11 


7.73 


AF e 


5.403 


0.431 


5.834 


0.074 


-1.098 


-0.709 


0.342 


1.569 


1.125 


2.169 


8.18 


6.82 , 6.88 


AG 


5.738 


0.737 


6.475 


0.114 


-0.891 


-0.465 


0.330 


1.351 


1.296 


2.945 


7.15 


6.98 



a Value derived for the q parameter (in units of cm s — ) obtained using the optimal calibration given bv. Kewlev &c Dopit al J2002J). 

Value derived for the q parameter (in units of cm s — ) obtained using the iterative procedure described in Kobulnicky &c Kcwlcv (2004). 
c For these models, for which wc derived 12+log(0/H)^8.4 following the T c method, we list the results of the empirical calibrations considering both the 
low and the high mctallicity branches. The <JkK04 parameters listed here for these models are for the low and high mctallicity branches, respectively. 



